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ABSTRACT 


A  high  pressure  apparatus  has  been  developed  for  measuring 
refractive  indices  and  compositions  of  coexisting  phases  in  vapor- 
liquid  equilibrium.  The  mixture  to  be  studied  is  contained  between  two 
pistons  which  are  connected  hydraulically  to  an  opposed  transmission 
motorized  Ruska  pump.  This  allows  the  pistons  to  be  moved  up  and  down 
simultaneously  while  maintaining  the  volume  between  them  constant.  In 
the  center  of  the  cell  a  pyrex  window  and  a  steel  mirror  form  a  30° 
prism  shaped  cavity.  The  refractive  index  of  the  fluid  in  the  cavity 
is  measured  by  knowing  the  prism  angle  and  measuring  the  deviation  angle 
of  a  beam  of  light  from  an  autocol 1 imati ng  telescope.  The  optical  system 
design  is  such  that  the  only  moving  part  is  the  telescope  and  the  devia¬ 
tion  angle  is  automatically  the  minimum  deviation  angle.  Small  samples 
are  taken  from  either  of  two  micrometering  sampling  valves  and  analyzed 
by  gas  chromatography . 

The  equipment  and  method  were  evaluated  by  comparison  of  the 
molar  volumes  calculated  from  the  refractive  index  with  those  obtained 
by  direct  volumetric  measurement  for  the  n-butane-carbon  dioxide  system 
at  100°F. 

The  compositions  and  refractive  indices  of  the  coexisting  vapor 
and  liquid  phases  of  the  i -butane-carbon  dioxide  and  i -butane-ethane 
systems  were  determined  in  the  temperature  range  100°F  to  250°F.  The 
pressure  range  extended  from  the  vapor  pressure  of  i -butane  to  the 
critical  regions  of  the  systems.  The  molar  volumes  of  the  phases  were 
calculated  from  their  measured  refractive  indices  and  their  molar  average 
refracti vi ties .  The  thermodynamic  consistency  of  the  data  was  tested. 


The  refractive  indices  of  carbon  dioxide,  ethane  and  liquid 
i -butane  were  determined  in  the  temperature  range  100°  to  250°F  up  to 
a  pressure  of  1500  psia. 
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I  -  INTRODUCTION 


Data  on  the  composition  of  equilibrium  vapor  and  liquid 
phases  of  hydrocarbon  and  many  other  systems  have  been  obtained  by 
several  different  methods  in  a  variety  of  different  experimental 
cells.  Only  rarely,  however,  have  data  on  phase  compositions  been 
obtained  simultaneously  with  density  data.  Aside  from  being  of  interest 
in  themselves,  the  coexisting  phase  density  data  are  useful  in  several 
other  ways.  They  are  a  valuable  asset  in  the  evaluation  of  vapor- 
liquid  equilibrium  prediction  methods.  With  equations  of  state,  for 
example,  the  phase  density  data  are  used  to  calculate  the  component 
fugacities  and  partial  molar  volumes.  Also,  the  phase  density  data 
are  an  intrinsic  part  of  the  enthalpy  and  entropy  departure  calculations. 

The  residual  transport  properties,  residual  tiiermal  conduct¬ 
ivity  and  residual  viscosity,  appear  to  be  well  correlated  to  density. 
Therefore,  knowledge  of  coexisting  phase  densities  can  be  used  to 
estimate  their  viscosities  and  thermal  conductivities.  This  work  was 
undertaken  therefore  to  produce  a  convenient  and  accurate  method  for 
obtaining  both  composition  and  density  data  at  elevated  pressures  on 
coexisting  vapors  and  liquids. 

The  refractive  index  was  chosen  as  a  means  for  accomplishing 
this  because  this  property  can  be  measured  easily  and  accurately,  and 
it  is  related  quantitatively  to  the  density.  Accordingly,  a  cell  was 
designed  which  was  capable  of  attaining  equilibrium  and  of  allowing 
the  compositions  and  the  refractive  indices  of  the  phases  to  be  measured 
with  convenience  and  accuracy. 

The  n-butane-carbon  dioxide  system  at  100°F  was  chosen  to 
evaluate  the  performance  of  the  cell  and  the  method.  Two  previously 


1 
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unreported  systems,  the  i -butane-ethane  system  and  the  i -butane-carbon 
dioxide  system,  were  then  studied  in  the  temperature  range  100°F  to 
250°F  at  pressures  extending  from  the  vapor  pressure  of  i -butane  to 
the  critical  regions  of  the  systems.  The  densities  of  the  phases  were 
calculated  from  their  measured  compositions  and  refractive  indices. 

The  cell  may  also  be  used  to  determine  single  phase  refractive 
indices.  The  refractive  indices  of  carbon  dioxide,  ethane  and  liquid 
i -butane  were  determined  in  the  temperature  range  100°  to  250°F  up  to  a 
pressure  of  1500  psia. 


II  -  RELATIONSHIP  OF  REFRACTIVE  INDEX  TO  DENSITY 


For  a  given  substance  the  refractive  index  is  a  property 

which  is  closely  related  to  the  density.  The  first  formula  expressing 

this  relationship  was  put  forward  by  Newton  in  1704,  and  since  then  a 

number  of  empirical  and  semi -empi ri cal  formulae  have  been  proposed. 

1  2 

Both  Batsanov  and  Partington  have  reviewed  this  historical  develop- 

3 

ment.  Hill,  Vaughan,  Price  and  Davies  have  presented  a  comprehensive 
treatment  of  how  dielectric  properties  and  refractive  index  are 
related  to  molecular  behavior. 

A  list  of  the  more  important  refractivity  formulae  is  given 

below : 


1)  Newton-Lapl ace  (1704,1806) 

1  2 

—  •  (n  -  1 )  =  constant 

2)  Beer-Gl ads tone-Dale  (1853,1863) 

—  •  (n  -  1 )  =  constant 

p 

3)  Lorentz-Lorenz  (1880) 

1  n2  -  1 

—  •  —7T— —  =  constant 
p  n2  +  2 

4)  Edwards  (1894) 


1  n  -  1 
p  n 


1  .01 


1  .02 


1  .03 


1 .04 


3 


constant 


Sim 


--  ■ 
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5)  Zecchini  (1895) 


1 

p 


2  n 

n  -  1 

3  x  9 

n  +2 


constant 


6)  Eykmann  (1895) 


1 

p 


n 


-  1 


n  +  0.4 


constant 


7)  Lichtenecker  (1926) 


-  •  (log  n)  =  constant 


8)  Macdonald  (1926) 


1 

p 


n2  -  1 

n 


constant 


1  .05 


1  .06 


1 .07 


1  .08 


The  simplest  formula,  Beer-Gl  adstone-Dale ,  holds  reasonably 

well  over  a  wide  density  range.  This  is  exemplified  in  Figure  1  which 

shows  a  plot  of  refractive  index  versus  density  for  the  data  of 

4  5  6 

Michels  et  al  on  ethylene  ,  carbon  dioxide  ,  and  nitrogen  at  pressures 
up  to  2000  atm.  The  Lorentz-Lorenz  ref racti vi ty ,  however,  is  the 
soundest  expression  relating  the  refractive  index  to  the  density 
because  it  is  the  only  formula  which  can  be  related  to  the  polariz¬ 
ability  of  the  molecule. 

If  an  electric  field,  F,  acts  on  a  molecule  the  average 
moment  induced  will  be  proportional  to  the  field  provided  that  the 
field  is  not  extraordi nari ly  large. 
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FIG;  1  RELATIONSHIP  OF  REFRACTIVE  INDEX  TO  DENSITY  FOR  ETHYLENE,  CARBON  DIOXIDE  AND  NITROGEN 
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rn  —  ot-p  F  1.09 

The  constant  of  proportionality  is  called  the  polarizability  of  the 
molecule  and  is  the  sum  of  the  electronic,  the  atomic,  and  the 
orientation  polarizabilities. 


a 


T 


=  ae  + 


a  +  a 

a  o 


1.10 


The  electronic  contribution  arises  from  the  electric  field 

causing  a  displacement  of  the  electrons  relative  to  the  nucleus  in 

each  atom;  the  atomic  contribution  from  a  displacement  of  the  atomic 

nuclei  relative  to  one  another;  and  the  orientation  contribution  from 

the  tendency  of  the  electric  field  to  align  molecules  with  permanent 

dipole  moments  parallel  to  itself. 

Each  of  the  three  types  of  polarizability  is  a  function  of 

the  frequency  with  which  the  applied  field  is  altered.  In  a  steady 

field  all  three  types  are  at  their  equilibrium  values.  When  the  field 

is  altered  with  a  frequency  in  the  range  of  10^  to  10^  cycles  per 

second  the  orientation  polarizability  fails  to  attain  its  equilibrium 

value  and  at  higher  frequencies  no  longer  contributes  to  the  total 

13  14 

pol ari zabi 1 i ty .  At  infra-red  frequencies,  10  to  10  cycles  per 

second,  the  atomic  polarizability  fails  to  attain  its  equilibrium 

value  leaving  only  the  electronic  polarizability  effective  in  the 

14  15 

visible  and  ultraviolet  regions  of  the  spectrum,  10  to  10  cycles 
per  second. 

The  Lorentz-Lorenz  formulation  is  given  by  Equation  1.11 


7 


M  n2  -  1 
p  n2  +  2 


4nN 

3e  ae 


constant 


1.11 


Experimentally  the  Lorentz-Lorenz  refractivity  is  not  exactly 

constant  but  increases  slightly  with  increasing  density  passing  through 

a  maximum  and  declining  thereafter.  This  behavior  is  exemplified  in 

Figure  2  which  shows  a  plot  of  versus  density  for  Michels  and 
.  5 

Hamers  data  on  carbon  dioxide  at  99.7°C. 

Several  theories  have  been  put  forward  to  account  for  these 
variations.  For  anisotropic  molecules  Raman  and  Krishnan^  deduced 
that  there  would  be  an  asymmetric  distribution  of  the  polarizable 
matter  surrounding  any  given  molecule  in  the  medium,  rather  than  a 
spherical  one  as  assumed  in  Lorentz-Forenz  formulation.  As  a  conse¬ 
quence,  the  anisotropic  nature  of  the  polarization  field  will  increase 

with  the  molar  density  causing  a  progressive  variation  of  R^ . 

8  9 

Kirkwood  and  Yvon  showed  that  even  for  spherical  molecules  the 
polarization  field  fluctuates  with  the  thermal  motion  of  the  mole¬ 
cules  and  produces  a  fluctuating  induced  moment  the  magnitude  of 
which  increases  with  increasing  density. 

Mazur  and  Mandel^and  Jansen  and  Mazur  ^  used  a  quantum 
mechanical  approach  and  found  that  the  polarizability  should  vary  with 

density,  increasing  when  attractive  forces  are  dominant  and  decreasing 

12 

when  repulsive  forces  are  dominant.  Buckingham  and  Pople  also  used 
a  quantum  mechanical  approach  and  proposed  that  R^  could  be  repre¬ 
sented  by  a  virial  expansion  in  molar  volume,  namely 

rll  =  rll  +  B/v  +  c/v2 


1.12 


6.70 
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FIG.  2  VARIATION  OF  THE  LORENTZ-LORENZ  REFRACTIVITY  WITH  DENSITY  FOR 
CARBON  DIOXIDE  AT  99.7°C 
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where  R^L  is  the  limit  of  the  molar  refractivity  at  zero  density,  and 
B  and  C  are  the  second  and  third  virial  refractivity  coefficients  res- 

n  \  R 

pectively.  A  plot  of  versus  1  / V  of  Michels  and  Hamers' 

data  for  carbon  dioxide  is  shown  in  Figure  3,  where  for  carbon 
dioxide  was  obtained  from  Stoll's^  low  pressure  data.  B  and  C  are 
obtained  from  the  intercept  and  slope  respectively. 

In  the  visible  region  of  the  spectrum  the  effect  of  fre¬ 
quency  on  R^  can  be  adequately  represented  by  making  only  R°^  a 
function  of  wavelength,  since,  within  experimental  error,  B  and  C 
are  independent  of  wavelength.  Figure  4  shows  the  effect  of  wave¬ 
length  on  R°^  for  CC^. 
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FIG.  3  DETERMINATION  OF  THE  SECOND  AND  THIRD  REFRACTIVITY  VIRIAL 
COEFFICIENTS  FOR  CARBON  DIOXIDE  AT  99.7°C 


6.80 
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FIG.  4  VARIATION  OF  THE  CARBON  DIOXIDE  LOW  DENSITY  REFRACTIVITY  WITH  WAVELENGTH 


Ill  -  METHODS  OF  MEASURING  REFRACTIVE  INDEX 


There  are  four  means  by  which  the  refractive  index  of  a 
fluid  can  be  measured: 

a.  Spectrometry 

b.  Interferometry 

c.  Critical  angle  refractometry 

d.  Immersion 

Spectrometers 

A  spectrometer  is  used  to  measure  the  minimum  angle  through 
which  a  beam  of  monochromatic  light  is  deviated  by  a  prism  shaped 
sample  of  the  fluid.  Figure  5  shows  a  prism  in  the  position  of  mini¬ 
mum  deviation.  In  this  position  the  angle  of  incidence  and  the  angle 
of  emergence  are  equal : 


i  .  =  i  =  i  3.01 

1  o 

The  refraction  angles  are  equal  to  each  other  and  one  half  the  prism 
apex  angle,  A: 


r.  =  r  =  A/2  3.02 

i  o 

The  angle  of  incidence  is  equal  to  one  half  the  apex  angle  plus  the 
minimum  deviation  angle,  D: 

i  =  1/2  (A  +  D)  3.03 

From  Snell's  law  of  refraction 
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FIG.  5  SCHEMATIC  REPRESENTATION  OF  A  SPECTROMETER  IN  THE  POSITION  OF  MINIMUM  DEVIATION  ANGLE 
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n-j  sin  i  =  sin  r  3.04 

where  n-|  is  the  refractive  index  of  air  and  the  refractive  index 
of  the  prism  material  ,  it  follows  that 


n 


2 


sin  1/2(A  +  D) 

sin  1/2(A) 


3.05 


Interfe rometers 

An  interferometer  or  etalon  consists  of  two  thin  plates 
separated  by  a  distance  d.  A  slitted  beam  of  monochromatic  light 
of  wavelength  A  is  sent  to  the  i nterferometer  and  the  interference 
rings  produced  are  observed  in  the  focal  plane  of  the  objective  lens 
of  a  telescope. 

The  measurement  of  the  refractive  index  is  usually  carried 
out  in  the  following  manner.  The  space  between  the  plates  is  first 
evacuated  and  the  interference  pattern  observed.  The  gas  whose 
refractive  index  is  to  be  measured  is  slowly  let  into  the  region 
containing  the  etalon.  The  number  of  interference  fringe  shifts,  k, 
are  counted.  The  refractive  index  of  the  gas  between  the  plates  can 
thus  be  obtained  from  the  formula 

n  -  1  =  -~  3.06 


Critical  Angle  Refractometers 

The  Abbe  and  Pul f rich  refractometers  are  two  types  of  commer¬ 
cial  refractometer  which  use  the  angle  at  which  total  internal  reflection 
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occurs  for  light  passing  from  glass  to  a  liquid  on  its  surface.  From 
Snell's  law  and  Figure  6  it  follows  that: 


n 


3 


n  sin  B 
n2  sin  tt/2 


r sin  B 


3.07 


where  n^  is  the  refractive  index  of  the  liquid  n^  the  refractive  index 
of  the  glass  and  B  the  critical  angle. 

Immersion  Methods 

The  immersion  method  is  based  on  the  disappearance  of  an 
interface  between  two  immiscible  fluids  or  a  fluid  and  a  solid  when 
their  refractive  indices  are  equal.  If  the  refractive  index  of  one 
of  the  phases  is  known  then  the  other  is  also  known. 
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FIG.  6  SCHEMATIC  REPRESENTATION  OF  THE  CRITICAL  ANGLE  OF  REFRACTION 


IV  -  HIGH  PRESSURE  REFRACTOMETERS 


Several  high  pressure  refractometers  have  been  built  using 

the  minimum  deviation  angle  principle.  Examples  of  these  are  described 

14  15  16 

by  Lyons  and  Poindexter  ,  Poulter  et  al  ,  Smith  ,  Schmidt  and 
17  18 

Thomas  ,  and  Teague  and  Pings  .  In  all  cases  a  separate  collimator 
and  telescope  were  used  and  no  means  of  mixing  the  cell  contents  were 
employed. 

Interferometers  have  been  used  frequently  to  measure  refrac¬ 
tive  indices  at  high  pressures.  Examples  of  these  are  described  by 

19  20  21  22  PI 

Siertsema  ,  Ayres  ,  Phillips  ,  Himstedt  and  Wertheimer  ,  Eisele  , 

Opladen^,  Michels  et  al^’^,  Belonogov  and  Gorbunkov^5, Diller 

27,28  .  c, .  .  .  .29 

,  and  Sliwinski 

The  critical  angle  and  immersion  techniques  have  had  limited 

application  at  elevated  pressures.  Bellingham  and  Stanley  Limited 

(London,  England)  have  a  pipe  refractometer  which  is  capable  of  operat- 

30 

ing  at  pressures  up  to  150  psi.  Francis  used  the  immersion  principle 
to  obtain  refractive  indices  of  many  liquified  gases  at  elevated  pres¬ 
sures.  He  determined  the  iso-optic  temperature  of  the  liquified  gas 
and  an  immiscible  fluid,  such  as  brine,  sealed  in  a  thick-wall  glass 
tube.  The  refractive  index  of  the  Immiscible  fluid  was  known  or  could 
be  measured  easily  and  therefore  the  refractive  index  of  the  liquified 
gas  at  the  iso-optic  temperature  was  known. 

Several  unique  high  pressure  refractometer  designs  have  been 
described  in  the  literature.  The  optical  system  for  measuring  the 

refractive  index  of  a  fluid  in  a  capillary  tube  was  described  by 

31 

Melville  and  Watson  .  A  similar  type  of  system,  the  capillary  pro- 

32 

jection  method,  was  later  used  by  Deryagin  et  al  and  Konstantinov 
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et  al  .A  refractometer  consisting  of  a  coiled  pyrex  rod,  a  constant 
intensity  light  source,  and  a  photoelectric  detector  was  described  by 
Shaw  and  Johnson^.  The  loss  of  intensity  of  the  light  shone  in  one 
end  of  the  rod  and  measured  at  the  other  was  found  to  be  directly  pro¬ 
portional  to  the  refractive  index  of  the  fluid  in  which  the  coiled  rod 
was  immersed.  A  high  pressure  refractometer  design  employing  an  internal 
rotating  glass  prism  and  no  other  moving  parts  has  been  described  by 
Vukalovich  et  al^. 

In  this  study,  a  high  pressure  refractometer  was  designed 
based  on  the  minimum  deviation  angle  principle.  The  optical  system 
differs  somewhat  from  previous  designs  in  that  the  only  moving  part 
is  an  autocollimating  telescope  and  the  measured  angle  is  automatically 
the  minimum  deviation  angle. 


V  -  PREVIOUS  WORK 


A  number  of  experimental  studies  have  been  made  to  determine 

the  relationship  between  density  and  refractive  index.  These  studies 

have  frequently  involved  a  lighter  gaseous  substance  which  is  confined 

and  then  compressed  so  that  the  refractive  index  can  be  obtained  over 

a  range  of  densities.  Measurements  of  this  kind  have  been  carried  out 

on  hydrogen  ,  deuterium  ,  helium  ,  argon  ,  xenon  , 

krypton  ,  nitrogen  *  ,  oxygen  ,  air  ,  carbon  dioxide 

19-21 ,24,30,35,38,45-47^  ethylene4, 25,48, 49>  methane37,49>  tetra_ 

fluoromethane  ,  water  ,  ammonia  ,  ethane  ,  propane 

*  *  ,  butane  ,  hydrogen  sulfide  *  ,  and  l -pentane  . 

The  refractive  index  and  density  of  compounds  which  are 
liquids  at  or  near  ambient  conditions  has  been  studied  extensively. 

r  «7  CO 

A  good  review  of  these  has  been  presented  by  Timmermans37’  . 
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VI  -  MIXTURES  AND  THE  ADDITIVITY  OF  THE 
LORENTZ-LORENZ  REFRACTI VITIES 

The  assumption  that  the  Lorentz-Lorenz  refracti vities  are 

additive  in  mixtures  can  be  justified  by  the  near  additivity  found 

in  a  large  number  of  mixtures  and  by  the  manner  in  which  atomic  and 

radical  refracti vities  can  be  summed  to  give  molecular  refractivi ties . 

59 

Smyth,  Engel,  and  Wilson  studied  the  refractive  index  and  density 

of  twelve  binary  liquid  mixtures.  Their  study  included  systems  such 

as  carbon  tetrachloride-acetone  and  n-heptane-ethyl  iodide,  covering 

their  full  compositional  ranges.  In  every  case  the  maximum  deviation 

from  additivity  was  less  than  0.1  percent.  Bloom  and  Rhodes6^*  showed 

that  the  refraction  of  molten  salts  are  close  to  the  refraction  of 

the  salts  in  aqueous  solution  and  also  that  the  refraction  of  molten 

mixtures  follow  the  law  of  additivity  within  experimental  error, 
fil 

Young  and  Finn  found  that  in  glasses  the  refracti vities  are  addi¬ 
tive  to  within  t  0.1  percent. 

62 

Recently  Burfield,  Richardson  and  Guereca  used  dielectric 

constant  measurements  and  a  molar  average  Clausius-Mosotti  function 

o 

(the  dielectric  constant  replaces  n  in  Equation  1.03)  to  calculate 
coexisting  phase  densities  in  the  helium-carbon  dioxide  system. 

Their  results  compared  favorably  with  previous  volumetric  measurements. 

Keilich  in  a  theoretical  study  of  molar  refracti vities  of 
dense  mixtures  has  shown  that  is  strictly  additive  only  in  the 
case  of  perfect  gases.  The  non-additive  terms  were  found  to  be  func¬ 
tions  of  the  polarizabilities,  molecular  symmetry  and  dipole  moment 
of  the  molecules  and  the  density  of  their  mixture. 
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Although  the  existence  of  non-additive  terms  appears  to  be 
well  founded  in  theory,  experimentally  these  terms  have  been  shown  to 
be  rather  small.  In  the  study  of  high  pressure  phase  equilibria  the 
compositions  of  the  phases  can  usually  be  measured  to  within  t  0.3 
mole  percent,  with  some  improvement  at  the  extreme  composition  ends. 
When  the  refractive  index  is  used  to  calculate  the  equilibrium  phase 
molar  volumes,  this  and  the  associated  errors  in  the  measurement  of 
temperature  and  pressure  will  make  the  error  due  to  the  non-additive 
terms  small  in  comparison.  Also, the  measurement  of  the  refractive 
indices  of  the  equilibrium  phases  can  be  made  conveniently  and 
accurately  which  is  generally  not  the  case  when  direct  volumetric 
measurements  are  attempted. 


VII  -  EQUIPMENT  DESIGN 


A  schematic  diagram  of  the  apparatus  and  its  associated 
equipment  is  shown  in  Figure  7  and  the  essential  features  of  the 
cell  design  are  shown  in  Figure  8.  It  has  an  overall  length  of  48 
cm.  and  was  machined  from  a  10  cm.  diameter  type  316  stainless  steel 
cylinder.  The  cell  proper  consists  of  three  parts,  two  cylinder- 
piston  end  sections  and  a  central  windowed  section.  The  three  sec¬ 
tions  are  bolted  together  with  the  high  pressure  seals  between  them 
made  by  flattened  teflon  0-rings.  Each  piston  has  a  10  cm.  travel 
and  is  confined  to  its  respective  cylinder.  The  function  of  the 
pistons  is  to  isolate  the  cell  contents  from  the  hydraulic  fluid, 
a  low  viscosity  silicone  oil,  and  to  provide  a  means  for  varying 
the  cell  volume.  The  piston  seal  is  effected  by  four  0-rings,  one 
teflon  and  three  neoprene. 

The  horizontal  section  A-A  through  the  center  of  the  cell 
shows  how  the  window,  a  pyrex  glass  disk  2.9  cm.  in  diameter  and  2.5 

,  _4 

cm.  thick  with  its  faces  flat  and  parallel  to  T  10  cm.,  and  a  type 
316  stainless  steel  mirror  form  the  boundaries  of  a  30°  prism.  The 
window  and  mirror  are  held  in  place  by  bolted  cover  plates  and  sealed 
against  lapped  surfaces  with  1.5  mm.  thick  glass-filled  teflon  gaskets. 
A  316  stainless  steel  sheathed  iron-constantan  thermocouple  has  its 
tip  exposed  to  the  fluid  forming  the  prism. 

A  second  horizontal  section  B-B  shows  one  of  the  sampling 
valves,  the  pressure  transducer  port  and  the  circulation  line  outlet. 
The  sampling  valves  are  internally  mounted  Circle  Seal  MV-92  micro¬ 
metering  valves  with  provision  made  for  flushing  out  and  evacuation 
of  the  low  pressure  sampling  lines.  The  circulation  system  used  to 
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SCHEMATIC  DIAGRAM  OF  THE  EQUILIBRIUM  CELL  AND 
ITS  ASSOCIATED  EQUIPMENT 
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FIG.  8  ESSENTIAL  FEATURES  OF  THE  EQUILIBRIUM  CELL  DESIGN 
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attain  equilibrium  consists  of  an  internal  check  valve  in  the  central 
section  and  an  external  3.3  mm.  diameter  stainless  steel  line  with  a 
shut-off  valve. 

Cel  1  Temperature  Control 

The  temperature  of  the  equilibrium  cell  is  maintained  by 
two  15  cm.  diameter  aluminum  environmental  control  shrouds  which  are 
placed  over  the  ends  of  the  cell  as  shown  in  Figure  8.  Each  shroud 
has  four  150  watt  strip  heaters  and  a  6  meter  cooling  coil.  A  thermo¬ 
couple  proportional  band  temperature  controller  (Thermoelectric  400) 
operates  the  heaters  in  both  shrouds  and  controls  the  cell  temperature 
to  within  -  .05°C  of  the  setpoint. 

Operational  Procedure 

The  two  pistons  are  hydraulically  driven  by  a  motorized  dual 
cylinder  high  pressure  pump  with  an  opposed  transmission  (Ruska  model 
2248  WII).  This  allows  both  pistons  to  be  move  simultaneously  up  or 
down  maintaining  the  volume  between  them  constant.  The  possible  work¬ 
ing  volume  ranges  from  10  cc.  to  175  cc. 

Mixing  and  equilibration  are  attained  by  running  the  pistons 
up  and  down.  On  the  upstroke,  the  internal  check  valve  closes  and  the 
fluids  are  forced  through  the  external  circulation  line  and  sprayed 
into  the  upper  cylinder.  On  the  downstroke  the  circulation  line  valve 
is  closed  and  the  fluids  pass  through  the  check  valve.  Five  to  ten 
cycles  are  usually  sufficient  to  attain  equilibrium. 

The  point  at  which  equilibrium  is  achieved  is  obvious  since 
the  reflected  light  image  becomes  stable  in  each  phase.  The  refractive 
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index  of  each  phase  is  easily  measured  by  moving  it  into  the  prism 
and  measuring  the  deviation  angle,  y.  The  temperature  is  measured 
with  an  i ron-constantan  316  stainless  steel  sheathed  thermocouple 
with  its  tip  exposed  to  the  fluid  in  the  prism  and  the  pressure  is 
measured  with  a  strain  gauge  pressure  transducer  (Consolidated 
Electronics)  calibrated  against  the  vapor  pressure  of  carbon  dioxide 
at  20°C,  and  is  believed  known  to  t  0.2  atm. 

Samples  of  liquid  and  vapor  are  expanded  to  a  pressure  of 
0.2  atm  through  the  micro-metering  valves  into  the  evacuated  line 
connected  to  the  gas  chromatograph  gas  sampling  valve.  The  pressure 
in  the  sampling  line  is  measured  with  a  differential  pressure  trans¬ 
ducer  (Pace  Wianko).  Two  samples  of  each  phase  are  taken  and  tripli¬ 
cate  chromatographs  are  run  on  each  sample.  The  average  size  of  each 

_3 

sample  is  10  g-moles  corresponding  to  a  0.2  percent  depletion  of  an 
average  load.  The  pressure  drop  associated  with  a  sample  taking  was 
always  less  than  0.1  atm.  The  gas  chromatograph  (Hewlett-Packard 
Model  700)  equipped  with  a  thermal  conductivity  cell  and  an  appropri¬ 
ate  column  was  calibrated  with  the  pure  components  over  a  sample  loop 
pressure  range  of  0.05  to  0.7  atm.  The  sample  pressures  were  adjusted 
so  that  the  peak  areas  of  the  components  were  within  the  calibrated 
intervals . 

The  hydraulic  fluid  is  preconditioned  before  it  enters  the 
cell  by  2.5  meters  of  high  pressure  tubing  which  has  been  coiled  into 
each  shroud.  The  whole  apparatus  is  insulated  by  a  double  layer  of 


27 


Armoflex  foam  rubber  insulation. 

Optical  System 

A  telescope  (Gaertner  M523  aperture  28  mm.,  f.l.  250  mm.) 

with  an  Abbe-Lamont  autocol limating  eyepiece  (Gaertner  L372)  is 

mounted  on  a  25  cm.  diameter  precision  rotary  table  (Karl  Kneise 

RTHP-10)  with  a  two  second  vernier  hand  wheel  and  a  maximum  table 

error  of  10  seconds  in  360  degrees  of  rotation.  The  table  is  centered 

under  a  point  on  the  front  surface  of  the  pyrex  cell  window.  Mono- 

o 

chromatic  light  at  a  wavelength  of  6328  A  is  provided  by  a  Helium- 
Neon  gas  laser  (Spectra  Physics  Model  132)  and  is  transmitted  to  the 
eyepiece  by  a  glass  fiber  light  guide. 

Figure  g  is  a  schematic  representation  of  the  angle  measure¬ 
ment  involved  in  the  determination  of  the  refractive  index.  The  angle 
Y  is  measured  by  autocol 1 imation  and  from  Snell's  law  it  follows  that 


sin 

Y 

=  n^  sin 

3 

7.01 

sin 

3 

=  n^  sin 

6 

7.02 

= 

"i 

sin  y/sin 

6 

7.03 

where  n-j  is  the  refractive  index  of  the  ambient  air,  n^  the  refractive 

index  of  the  glass,  and  n ^  the  refractive  index  of  the  fluid  forming 

prism.  The  angle  6  is  the  prism  angle  and  is  equal  to  the  measured 

angle  y  when  n-j  =  n^,  that  is,  the  cell  contains  ambient  air.  The 

measured  angles  have  a  repeatability  of  -  4  seconds  and  a  table  error 

of  -  2  seconds  over  the  measured  range.  The  refractometer  has  an  over- 

+  -5 

all  possible  error  of  -  6  x  10  in  the  refractive  index  and  a  possible 
range  of  1.0  to  1.7. 
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This  optical  system  has  three  advantages  over  a  conventional 
spectrometer  which  normally  consists  of  a  60°  prism  with  a  separate 
collimator  and  telescope^. 

1.  The  time  required  to  make  a  measurement  is  reduced  from  the  usual 
five  to  ten  minutes  to  less  than  one  minute  because  the  minimum 
deviation  angle  is  obtained  without  trial  and  error  and  without 
movement  of  the  prism,  in  this  case  the  cell. 

2.  Only  one  high  pressure  window  with  flat  and  parallel  faces  is 
required  as  opposed  to  two  of  equal  thickness. 

3.  The  rotary  table  can  be  centered  under  any  point  on  the  front  sur¬ 
face  of  the  glass  window  rather  than  the  optical  center  of  the  prism 
as  in  a  conventional  spectrometer. 


VIII  -  EQUIPMENT  EVALUATION 


The  carbon  dioxide-n-butane  binary  was  studied  at  100. 0°F 

65 

and  compared  with  the  data  of  Olds,  Reamer,  Sage  and  Lacey  and 

fifi 

Poettmann  and  Katz  .  The  n-butane  used  in  this  study  was  Matheson 

instrument  grade.  A  chromatographic  analysis  of  this  material  showed 

it  to  contain  99.5  mole  percent  n-butane,  0.4  mole  percent  propane, 

0.03  mole  percent  i -butane  and  0.07  mole  percent  nitrogen  and/or 

methane.  The  carbon  dioxide  was  obtained  from  Canadian  Liquid  Air 

Ltd.  and  contained  better  than  99.9  mole  percent  carbon  dioxide. 

Figure  10  shows  the  experimental  compositional  isothermal 

data  and  those  of  the  above  authors.  It  will  be  noted  that  the  present 

65 

data  and  the  data  of  Olds  et  al  show  some  significant  differences  when 

66 

compared  to  the  earlier  work  of  Poettmann  and  Katz  . 

o 

Figure  11  shows  the  refractive  index  (6328A)  of  n-butane, 
carbon  dioxide,  and  the  coexisting  liquid  and  gas  phases  of  their  mix¬ 
ture  at  100. 0°F  in  the  pressure  range  4  to  80  atm.  In  order  to  calcu¬ 
late  the  density  from  the  refractive  index  and  the  composition  the 
molar  ref racti vi ties  of  the  pure  components  were  required. 

The  refractive  index  of  carbon  dioxide  has  been  studied  by 

5 

Michels  and  Hamers  in  the  temperature  range  25  to  100°C  and  at  pres- 

21 

sures  up  to  2400  atm.,  and  by  Phillips  at  34°C  at  pressures  up  to 

1 3 

100  atm.  Stoll  determined  the  refractive  index  of  carbon  dioxide 
at  20°C  and  1  atm.  pressure  for  wavelengths  from  the  ultraviolet  to 
the  infrared.  The  first  term  in  the  vi ri al  expansion  of  Equation  1.12 
for  Rll  of  carbon  dioxide  at  6328ft  was  obtained  from  Stoll's  data  and 
the  second  and  third  terms  were  obtained  from  a  plot  of  CRll~Rll° ) V 
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FIG.  10  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION  DIAGRAM  FOR  THE  n-BUTANE-CARBON 
DIOXIDE  SYSTEM  AT  100°F 
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PRESSURE-REFRACTIVE  INDEX  DIAGRAM  FOR  THE  n-BUTANE-CARBON  DIOXIDE  SYSTEM 
AT  100°F 
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versus  1/V  of  Michels  and  Hamers'  high  density  data.  The  lower  density 
data  could  not  be  used  because  of  the  excessive  scatter  on  this  very 
sensitive  plot. 

For  n-butane,  several  measurements  have  been  reported,  all  at 
°  53 

5893A.  Grosse  determined  the  refractive  index  of  the  liquid  at  low 

temperatures  and  atmospheric  pressure  using  a  modified  Abbe  refracto- 

30 

meter.  Francis  determined  the  iso-optic  temperature  of  the  liquified 

gas  and  an  immiscible  fluid  or  solution  for  which  the  refractive  index 

29 

was  known  as  a  function  of  temperature.  Sliwinski  measured  the 
refractive  index  and  density  of  the  coexisting  liquid  and  vapor  phases 
in  the  temperature  range  0°  to  94°C. 

The  refractive  index  of  carbon  dioxide  and  n-butane  were 
measured  at  100. 0°F  and  pressures  up  to  100  atm.  The  molar  refractivity 
of  n-butane  was  calculated  using  molar  volumes  obtained  from  Sage 
Webster  and  Lacey6^.  The  smoothed  data  are  given  in  Table  1. 

The  mixture  refractivity,  (R|_L)m  was  calculated  by  using  a 
molar  average  mixing  rule 


“  lVRLL>i  8'01 

where  x^  is  the  mole  fraction  of  each  component.  The  molar  volume  of 
the  mixture  can  then  be  obtained  from  the  refractive  index  of  the  phase 
and  the  mixture  refractivity.  The  experimental  refractive  indices 

o 

(6328A) ,  the  compositions,  and  the  calculated  refracti vities  and  molar 
densities  of  the  coexisting  phases  in  the  carbon  dioxide-n-butane 
binary  are  given  in  Table  2. 

Figure  12  shows  the  calculated  densities  of  the  coexisting 
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TABLE  1 

REFRACTIVE  INDICES  OF  CARBON  DIOXIDE  AND  n-BUTANE  AT  100. 0°F 


n-Butane  Carbon  Dioxide 


Pressure 

Refractive 

LL 

Refractive 

LL 

atm. 

index 

cc/q-mole 

index 

cc/q-mole 

10 

1.3195 

20.44 

1 .0042 

6.64 

20 

1.3205 

20.43 

1 .0086 

6.65 

30 

1.3215 

20.42 

1.0136 

6.65 

40 

1.3226 

20.41 

1 .0192 

6.65 

50 

1.3236 

20.40 

1.0265 

6.65 

60 

1.3247 

20.39 

1 .0353 

6.64 

70 

1.3258 

2Q.38 

1.0478 

6.64 

80 

1.3267 

20.37 

1 .0738 

6.64 

90 

1.3277 

20.36 

1.1393 

6.61 

100 

1.3287 

20.35 

1 .1565 

6.60 

o 

a  relative  to  vacuum  at  6328A. 


Table  2.  Coexisting  Vapor  and  Liquid  Phase  Properties  in  the  Carbon  Dioxide-n-Butane  System  at  100. 0°F 
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FIG.  12  PRESSURE-EQUILIBRIUM  PHASE  DENSITY  DIAGRAM  FOR  THE  n-BUTANE-CARBON  DIOXIDE 
SYSTEM  AT  100°F 
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liquid  and  vapor  phases  and  also  the  data  of  Olds  et  al65  which  were 
obtained  by  direct  volumetric  measurement  for  the  binary  system  at 
100°F .  The  average  deviation  of  the  calculated  density  minus  that 
of  Olds  is  +1.1  percent  for  the  liquid  and  -1.8  percent  for  the  vapor. 

The  standard  deviation  of  all  the  points  is  t  1.5  percent. 

The  molar  refractivi ties  of  the  liquid  and  vapor  phases  were 
calculated  from  the  compositions  of  the  phases  using  the  ref racti vi ties 
of  the  component  substances  in  their  natural  state  at  the  operating 
temperature  and  pressure.  For  example,  at  100°F  and  40  atm.  the  carbon 
dioxide  exists  as  a  gas  and  has  a  refractivi ty  of  6.65  cc/g-mole  while 
the  n-butane  exists  as  a  liquid  and  has  a  refractivi ty  of  20.41  cc/g-mole. 
The  molar  ref racti vi ties  of  the  coexisting  liquid  and  vapor  were  cal¬ 
culated  using  these  values.  Having  obtained  the  refracti vi ties  sub¬ 
stitution  of  the  refractive  indices  yields  their  molar  densities. 

In  the  section  which  follows  the  same  procedure  is  used  to 
calculate  the  molar  densities  of  the  coexisting  phases  for  the  i -butane- 
carbon  dioxide  system  and  the  i -butane-ethane  system.  The  pure  component 
data  can  be  obtained  from  the  experimental  data  listed  in  Appendix  D  or 
can  De  calculated  using  the  refractivity  virial  coefficients  listed  in 
Appendix  E. 


IX  -  EXPERIMENTAL  DATA 


Two  binaries  are  reported  in  this  work,  the  i -butane-ethane 
system  and  the  i -butane-carbon  dioxide  system.  The  i -butane-ethane 
system  has  recently  been  studied  by  Skripa  et  al  in  the  temperature 
range  -94° F  to  32° F.  The  i -butane-carbon  dioxide  system  has  not 
previously  been  reported  in  the  literature. 

In  this  work  both  of  the  above  systems  were  studied  at  100°, 
160°,  220°,  and  250°F  at  pressures  from  the  vapor  pressure  of  i -butane 
to  the  critical  region  of  the  system.  The  refractive  indices  of  the 
coexisting  phases  were  measured  and  the  molar  volumes  calculated  using 
a  molar  average  Lorentz-Lorenz  refractivity.  The  refractive  indices 
of  i -butane,  ethane,  and  carbon  dioxide  were  measured  at  100°,  160°, 
220°  and  250°F  at  pressures  from  atmospheric  to  1500  psia. 

Apparatus 

The  apparatus  used  was  the  windowed  double  piston  cell  des¬ 
cribed  in  Chapter  VII.  It  was  found  that  the  hydraulic  fluid  preheating 
coils  were  not  long  enough  for  operations  at  temperatures  above  200°F. 
The  oil  was  therefore  initially  passed  through  a  coil  immersed  in  a 
200°F  bath  before  entering  the  preheat  coils  in  the  shroud.  The  large 
temperature  difference  between  the  oil  in  the  pump  and  the  oil  in  the 
cell  caused  a  small  change  in  the  working  volume  of  the  cell  when  the 
pistons  were  being  moved  simultaneously  up  or  down.  This  situation 
was  remedied  by  connecting  a  10  cc.  displacement  pump  to  one  of  the 
oil  lines.  The  pressure  in  the  cell  was  easily  maintained  constant 
manual ly  to  0.2  psi . 
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The  temperature  was  measured  with  an  i ron-constantan  thermo¬ 
couple  sheathed  in  316  stainless  steel  with  its  reference  junction  in 
an  ice  bath  and  its  measuring  junction  in  the  cell  contents.  The 
thermocouple  was  calibrated  at  the  triple  point  and  steam  point  of 
water.  The  cell  temperature  was  controlled  to  within  t  0.1  °F  of  the 
set  operating  temperature. 

The  pressure  was  measured  with  a  0-1500  psi  pressure  trans¬ 
ducer  calibrated  at  atmospheric  pressure  and  the  vapor  pressure  of 
carbon  dioxide  at  20°C.  The  maximum  combined  non-linearity  and 
hysteresis  of  the  transducer  was  t  3  psi. 

The  phase  compositions  were  determined  by  gas  chromatography 
as  described  earlier.  The  thermal  conductivity  cell  was  maintained  at 
200°C  and  150  ma  current  with  a  Helium  flow  rate  of  25  cc/min.  For 
the  i -butane-carbon  dioxide  system  a  six  foot  10  percent  UC-W98  80-100S 
column  maintained  at  20°C  was  used.  For  the  i -butane-ethane  system  a 
six  foot  Porapak  Q  column  maintained  at  100°C  was  used. 

The  refractive  indices  of  the  phases  were  determined  by  measur¬ 
ing  the  minimum  deviation  angle  of  a  prism  of  the  fluid  as  described 
earl ier. 

Materi als 

Matheson  instrument  grade  i-butane  was  used.  It  was  analyzed 
and  found  to  contain  99.9  mole  %  i-butane,  the  major  impurity  being 
propane.  The  ethane  used  was  Phillips  research  grade  with  a  reported 
analysis  of  99.94  mole  %  ethane.  The  carbon  dioxide  was  obtained  from 
Canada  Liquid  Air  Ltd.  It  was  analyzed  and  found  to  contain  better  than 
99.9  mole  %  carbon  dioxide. 


40 


Resul ts 

A.  Refractive  Indices  of  Ethane,  Carbon  Dioxide  and  i -Butane 

The  smoothed  refractive  indices  of  ethane,  carbon  dioxide, 

and  i -butane  at  pressures  up  to  1500  psia  in  the  temperature  range 

100°F  to  250°F  are  given  in  Tables  3,  4  and  5,  respectively.  Figures 

13,  14,  and  15  show  the  experimentally  determined  refractive  index 

data,  which  are  also  tabulated  in  the  Appendix.  The  molar  refracti vi ties 

of  i -butane  were  calculated  by  using  the  volumetric  data  of  Sage  and  Lacey^9 

and  the  measured  refractive  indices.  The  molar  refracti vi ties  used  for 

29  0 

ethane  were  those  reported  by  Sliwinski  with  R^  corrected  to  the 
He-Ne  lazer  frequency.  For  carbon  dioxide  the  values  obtained  from 

5 

Michels  and  Hamers  high  density  data  were  used. 

B.  Coexisting  Phase  Properties  of  the  i-Butane-Carbon  Dioxide  System 

Plots  of  the  experimentally  determined  phase  compositions  and 
refractive  indices  and  the  calculated  volumetric  data  for  the  i -butane- 
carbon  dioxide  system  are  shown  in  Figures  16,  17  and  18,  respectively. 

These  data  are  also  tabulated  in  the  Appendix.  The  smoothed  data  are 
given  in  Table  6  and  the  K-factor  versus  pressure  curves  are  shown  in 
Figure  19. 

C.  Coexisting  Phase  Properties  of  the  i -Butane-Ethane  System 

Plots  of  the  experimentally  determined  phase  composition, 
refractive  index  and  volumetric  data  for  the  i -butane-ethane  systems 
are  shown  in  Figures  20,  21  and  22,  respectively.  These  data  are  also 
tabulated  in  the  Appendix.  The  smoothed  data  are  given  in  Table  7  and 
the  K-factor  versus  pressure  curves  are  shown  in  Figure  23. 
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REFRACTIVE  INDEX,  6328A 

FIG.  13  EFFECT  OF  TEMPERATURE  AND  PRESSURE  ON  THE 
REFRACTIVE  INDEX  OF  ETHANE 
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Table  3.  Effect  of  Temperature  and  Pressure  on  the 
Refractive  Index  of  Ethane 


Pressure  Refractive  Index* 


psia 

100. 3°F 

160. 4°F 

219.6 

249.6 

100 

1 .0047 

1.0043 

1.0039 

1.0037 

200 

1.0100 

1 .0089 

1.0080 

1.0075 

300 

1.0159 

1.0138 

1.0122 

1.0115 

400 

1 .0230 

1.0191 

1.0166 

1.0157 

500 

1.0315 

1.0251 

1.0212 

1.0199 

600 

1.0438 

1 .0318 

1.0264 

1.0245 

700 

1.0611 

1.0395 

1.0319 

1.0295 

750 

1.0788 

1.0440 

1.0348 

1.0321 

800 

1 .1365 

1.0490 

1.0378 

1.0348 

850 

1.1620 

1.0540 

1.0410 

1.0375 

900 

1.1712 

1.0591 

1.0443 

1.0402 

1000 

1 .1835 

1.0716 

1.0514 

1.0461 

1100 

1.1916 

1.0865 

1.0587 

1.0521 

1200 

1.1979 

1.1023 

1.0665 

1.0585 

1300 

1.2025 

1.1180 

1.0747 

1.0652 

1400 

1 .2063 

1.1321 

1.0831 

1.0723 

1500 

1.2100 

1.1440 

1.0919 

1.0795 

o 

*Relative  to  vacuum  at  6328A 
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1400 
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1.16 


REFRACTIVE  INDEX,  6328A 

FIG.  14  EFFECT  OF  TEMPERATURE  AND  PRESSURE  ON  THE  REFRACTIVE 
INDEX  OF  CARBON  DIOXIDE 
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Table  4.  Effect  of  Temperature  and  Pressure  on  the 
Refractive  Index  of  Carbon  Dioxide 


Pressure 

psia 

100. 2°F 

Refractive 
160. 0°F 

Index* 

219. 8°F 

250. OF 

100 

1 .0027 

1.0025 

1.0023 

1.0021 

200 

1.0055 

1.0050 

1.0046 

1.0042 

300 

1 .0086 

1.0077 

1.0069 

1.0064 

400 

1.0120 

1.0104 

1 .0093 

1.0087 

500 

1.0158 

1.0134 

1.0118 

1.0110 

600 

1.0199 

1.0165 

1.0144 

1.0135 

700 

1 .0246 

1.0199 

1.0172 

1.0161 

800 

1.0299 

1.0235 

1.0201 

1.0187 

900 

1.0364 

1.0274 

1.0231 

1.0214 

1000 

1.0454 

1.0314 

1.0261 

1.0241 

1050 

1.0508 

1.0336 

1.0276 

1.0255 

1100 

1.0570 

1.0358 

1.0292 

1.0269 

1150 

1.0657 

1.0382 

1.0308 

1.0283 

1200 

1.0805 

1.0408 

1.0324 

1.0298 

1250 

1.1105 

1.0434 

1.0341 

1.0312 

1300 

1 .1333 

1.0462 

1.0358 

1.0327 

1350 

1.1438 

1.0492 

1.0376 

1.0342 

1400 

1.1499 

1.0520 

1.0395 

1.0358 

1500 

1  .1580 

1.0587 

1.0432 

1.0390 

o 

*Relative  to  vacuum  at  6328A 
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FIG.  15  EFFECT  OF  TEMPERATURE  AND  PRESSURE  ON  THE  REFRACTIVE 
INDEX  OF  LIQUID  I  -  BUTANE 
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Table  5.  Effect  of  Temperature  and  Pressure  on  the 
Refractive  Index  of  Liquid  i -Butane 


Pressure 

psia 

100. OF 

Refractive 
159. 8F 

Index* 

220. OF 

250. OF 

100 

1.3068 

- 

- 

- 

200 

1.3078 

1.2779 

- 

- 

300 

1 . 3088 

1 .2800 

- 

- 

400 

1.3097 

1 .2820 

1.2406 

- 

500 

1 .3107 

1.2838 

1.2449 

1  .2166 

600 

1 .3117 

1.2852 

1.2486 

1.2239 

700 

1 .3125 

1.2866 

1.2519 

1.2300 

800 

1 .3132 

1.2878 

1 .2548 

1.2349 

900 

1.3141 

1 . 2889 

1.2576 

1.2389 

1000 

1.3149 

1 .2901 

1.2599 

1.2422 

1100 

1 .3158 

1 .2913 

1.2620 

1.2452 

1200 

1 .3165 

1.2925 

1.2641 

1 .2480 

1300 

1.3173 

1.2937 

1.2661 

1 .2506 

1400 

1.3181 

1.2948 

1.2680 

1.2530 

1500 

1 .3190 

1.2960 

1.2700 

1.2553 

*Relative 

to  vacuum 

0 

at  6328A 

PRESSURE,  psia 
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MOLE  FRACTION  C02 

FIG.  16  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION 

DIAGRAM  FOR  THE  i-BUTANE-CARBON  DIOXIDE  SYSTEM 
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REFRACTIVE  INDEX,  6328A 


FIG.  17  PRESSURE-EQUILIBRIUM  PHASE  REFRACTIVE  INDEX  DIAGRAM 
FOR  THE  i-BUTANE-CARBON  DIOXIDE  SYSTEM 


PRESSURE,  psia 
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MOLAR  VOLUME,  cu  ft/lb  mole 


FIG.  18  PRESSURE-EQUILIBRIUM  PHASE  MOLAR  VOLUME 

DIAGRAM  FOR  THE  i-BUTANE-CARBON  DIOXIDE  SYSTEM. 
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Table  6.  Smoothed  Vapor-Liquid  Equilibria  Data  for  the 
i -Butane-Carbon  Dioxide  System 
100. 0°F 


Pressure 

Compos i tion 

Refractive 

Index 

Molar  Vol ume 

psia 

mole  frac 

.  co2 

cu. ft. /lb-mole 

Li  quid 

Vapor 

Liquid 

Vapor 

Liqui  d 

Vapor 

100 

0.021 

0.228 

1.3077 

1.0078 

1.74 

49.5 

150 

0.064 

0.477 

1.3031 

1.0094 

1.71 

34.0 

200 

0.107 

0.600 

1 .2982 

1.0112 

1.67 

26.0 

250 

0.150 

0.677 

1.2935 

1.0129 

1.64 

20.7 

300 

0.195 

0.723 

1 .2881 

1.0149 

1.60 

17.0 

350 

0.241 

0.758 

1.2829 

1.0169 

1.56 

14.1 

400 

0.287 

0.786 

1.2770 

1.0190 

1.53 

12.0 

500 

0.382 

0.821 

1.2644 

1 .0238 

1 .46 

9.10 

600 

0.482 

0.843 

1.2500 

1.0293 

1.40 

7.18 

700 

0.586 

0.862 

1.2331 

1 .0358 

1.35 

5.77 

800 

0.691 

0.879 

1.2138 

1.0437 

1.30 

4.68 

900 

0.782 

0.895 

1.1903 

1.0539 

1.28 

3.67 

1000 

0.854 

0.910 

1.1610 

1.0700 

1 .33 

2.70 

1050 

0.890 

0.916 

1.1409 

1.0935 

1 .44 

2.14 

1080 

0.918 

0.918 

1.110 

1.110 

1.68 

1.68 

160. 0°F 

Pressure 

Compos i tion 

Refracti  ve 

Index 

Mol  ar 

Vol  ume 

psi  a 

mole  frac 

.  co2 

cu. ft. /lb-mole 

Liquid 

Vapor 

Liquid 

Vapor 

Liquid 

Vapor 

200 

0.023 

0.154 

1 .2752 

1.0169 

1.90 

26.0 

250 

0.055 

0.303 

1.2719 

1.0189 

1.87 

21.1 

300 

0.088 

0.400 

1.2683 

1.0210 

1.85 

17.5 

350 

0.119 

0.471 

1.2643 

1.0231 

1.82 

14.7 

400 

0.149 

0.530 

1 .2603 

1 .0254 

1.80 

12.6 

450 

0.181 

0.569 

1.2562 

1.0279 

1.78 

11.0 

500 

0.213 

0.600 

1.2570 

1.0302 

1.76 

9.71 

600 

0.275 

0.645 

1 .2423 

1 .0357 

1.74 

7.77 

700 

0.338 

0.678 

1 .2320 

1.0420 

1.72 

6.38 

800 

0.403 

0.703 

1.2195 

1.0500 

1.72 

5.21 

900 

0.472 

0.718 

1.2033 

1.0606 

1.72 

4.24 

950 

0.513 

0.720 

1.1930 

1.0675 

1 .73 

3.78 

1000 

0.561 

0.720 

1 .1794 

1.0818 

1.76 

3.31 

1050 

0.618 

0.711 

1.1590 

1.0930 

1.88 

2.83 

1082 

0.672 

0.672 

1.129 

1.129 

2.23 

2.23 
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Table  6.  (continued) 
220°F 


Pressure 

Composition 

Refractive 

Index 

Molar 

Vol ume 

psia 

mole  frac 

.  co2 

cu. ft. /lb-mole 

Liquid 

Vapor 

Liquid 

Vapor 

Liquid 

Vapor 

350 

0.021 

0.084 

1.2373 

1 .0347 

2.20 

13.1 

400 

0.050 

0.179 

1.2340 

1.0372 

2.17 

11.5 

450 

0.078 

0.246 

1.2304 

1.0404 

2.13 

10.3 

500 

0.106 

0.307 

1 .2270 

1.0438 

2.11 

9.10 

550 

0.131 

0.349 

1.2235 

1.0475 

2.10 

8.18 

600 

0.159 

0.377 

1.2192 

1 .0516 

2.10 

7.33 

650 

0.187 

0.398 

1.2143 

1.0561 

2.11 

6.60 

700 

0.213 

0.410 

1.2095 

1.0611 

2.12 

5.98 

750 

0.240 

0.415 

1.2027 

1.0671 

2.15 

5.40 

800 

0.268 

0.416 

1.1947 

1.0742 

2.19 

4.87 

850 

0.300 

0.411 

1.1840 

1.0840 

2.26 

4.30 

900 

0.345 

0.403 

1 .1640 

1.1030 

2.42 

3.66 

925 

0.380 

0.380 

1.133 

1.133 

2.91 

2.91 

Pressure 
psi  a 

Composition 
mole  frac.  C02 

Liquid  Vapor 

250. 0°F 

Ref racti ve 

Liquid 

Index 

Vapor 

Molar  Volume 
cu. ft. /lb-mole 

Liquid  Vapor 

450 

0.015 

0.039 

1.2090 

1.0575 

2.51 

8.48 

500 

0.044 

0.113 

1 .2039 

1 .0618 

2.51 

7.69 

550 

0.072 

0.168 

1.1979 

1.0664 

2.52 

6.89 

600 

0.100 

0.202 

1.1908 

1.0722 

2.53 

6.04 

650 

0.128 

0.230 

1.1817 

1.0819 

2.59 

5.21 

700 

0.158 

0.237 

1 .1658 

1.0958 

2.77 

4.34 

738 

0.202 

0.202 

1 .129 

1.129 

3.36 

3.36 

EQUILIBRIUM  RATIO, 
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PRESSURE,  psia 

FIG.  19  EQUILIBRIUM  RATIOS  TOR  i-BUTANE  AND  CARBON  DIOXIDE 
IN  THE  i-BUTANE-CARBON  DIOXIDE  SYSTEM 
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FIG.  20  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION  DIAGRAM 
FOR  THE  i-BUTANE-ETHANE  SYSTEM 
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REFRACTIVE  INDEX,  6328A 


FIG.  21  PRESSURE-EQUILIBRIUM  PHASE  REFRACTIVE  INDEX 
DIAGRAM  FOR  THE  i- BUTANE-ETHANE  SYSTEM 
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MOLAR  VOLUME,  cu  ft/lb  mole 

FIG.  22  PRESSURE-EQULI BRIUM  PHASE  MOLAR  VOLUME  DIAGRAM 
FOR  THE  i-BUTANE-ETHANE  SYSTEM 
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Table  7.  Smoothed  Vapor  Liquid  Equilibria  Data 
for  the  i -Butane-Ethane  System 
100. 6°F 


Pressure 

Composi tion 

Refractive 

Index 

Molar 

Vol ume 

psia 

mole  frac. 

ethane 

cu. ft. /lb-mole 

Liquid 

Vapor 

Liquid 

Vapor 

Liquid 

Vapor 

100 

0.057 

0.334 

1.3045 

1.0083 

1 .71 

51.1 

150 

0.161 

0.545 

1.2970 

1.0111 

1.65 

33.9 

200 

0.261 

0.676 

1.2892 

1.0141 

1.61 

24.8 

250 

0.355 

0.755 

1.2811 

1.0171 

1.57 

19.2 

300 

0.443 

0.809 

1.2728 

1.0205 

1.54 

15.4 

350 

0.523 

0.840 

1.2640 

1.0243 

1.51 

12.7 

400 

0.597 

0.864 

1.2546 

1.0286 

1  .49 

10.6 

450 

0.667 

0.887 

1.2450 

1.0331 

1.48 

9.00 

500 

0.734 

0.907 

1.2343 

1 .0385 

1 .47 

7.57 

550 

0.794 

0.923 

1.2222 

1.0447 

1.49 

6.39 

600 

0.850 

0.940 

1.2083 

1.0522 

1 .53 

5.39 

650 

0.900 

0.952 

1.1915 

1 .0624 

1.60 

4.48 

700 

0.943 

0.971 

1 .1695 

1.0791 

1.74 

3.55 

740 

0.980 

0.980 

1.128 

1.128 

2.33 

2.33 

160. 4°F 

Pressure 

Composi tion 

Refractive 

Index 

Molar 

Vol ume 

psia 

mole  frac. 

ethane 

cu.  ft. /lb-moli 

Liquid 

Vapor 

Liquid 

Vapor 

Liqui  d 

Vapor 

200 

0.056 

0.213 

1 .2730 

1.0179 

1.89 

25.0 

250 

0.123 

0.370 

1.2678 

1.0211 

1.86 

19.6 

300 

0.189 

0.475 

1.2620 

1.0247 

1.84 

15.9 

350 

0.249 

0.550 

1.2555 

1 .0281 

1.82 

13.3 

400 

0.308 

0.603 

1 .2490 

1.0321 

1  .81 

11.3 

450 

0.364 

0.650 

1.2420 

1 .0364 

1.80 

9.75 

500 

0.420 

0.687 

1 .2341 

1.0411 

1.80 

8.51 

550 

0.475 

0.713 

1 .2260 

1.0465 

1.80 

7.45 

600 

0.530 

0.734 

1 .2165 

1.0526 

1 .81 

6.45 

650 

0.581 

0.751 

1.2063 

1.0599 

1.85 

5.54 

700 

0.633 

0.763 

1  .1942 

1.0690 

1.92 

4.68 

750 

0.683 

0.773 

1.1773 

1.0818 

2.02 

3.89 

800 

0.738 

0.778 

1 .1450 

1.1050 

3.05 

810 

0.760 

0.760 

1.126 

1.126 

2.60 

2.60 
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Table  7.  (continued) 


219. 7°F 


Pressure 

psia 

Composition 
mole  frac.  ethane 
Liquid  Vapor 

Refractive  Index 

Liquid  Vapor 

Molar  Volume 
cu. ft. /lb-mole 
Liquid  Vapor 

350 

0.037 

0.088 

1.2341 

1.0371 

2.20 

12.8 

400 

0.091 

0.197 

1.2284 

1.0416 

2.19 

10.8 

450 

0.144 

0.280 

1.2225 

1.0465 

2.19 

9.35 

500 

0.195 

0.348 

1.2160 

1.0520 

2.20 

8.16 

550 

0.242 

0.399 

1.2085 

1.0585 

2.22 

7.07 

600 

0.290 

0.433 

1.2000 

1.0663 

2.27 

6.08 

650 

0.334 

0.457 

1.1890 

1.0764 

2.32 

5.16 

700 

0.386 

0.466 

1.1700 

1.0960 

2.49 

4.10 

730 

0.450 

0.450 

1.13 

1.13 

3.14 

3.14 

Pressure 

psia 

Composi tion 
mole  frac.  ethane 
Liquid  Vapor 

249. 6°F 

Refracti  ve 

Liquid 

Index 

Vapor 

Molar  Volume 
cu. ft. /lb-mole 
Liquid  Vapor 

450 

0.013 

0.037 

1 .2090 

1.0570 

2.47 

8.80 

500 

0.064 

0.128 

1 .2045 

1.0645 

2.52 

7.42 

550 

0.115 

0.187 

1.1940 

1.0751 

2.64 

6.00 

600 

0.164 

0.217 

1.1725 

1.0950 

2.92 

4.70 

625 

0.200 

0.200 

1.13 

1.13 

3.52 

3.52 

EQUILIBRIUM  RATIO,  y/x 
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FIG.  23  EQUILIBRIUM  RATIOS  FOR  I  -  BUTANE  AND  ETHANE 
IN  THE  ETHANE-! -BUTANE  SYSTEM. 


X  -  THERMODYNAMIC  CONSISTENCY  AND  PREDICTION 


OF  THE  DATA 


Chueh-Prausni tz  Parameters 

The  smoothed  experimental  pressure-composition  data  were 
processed  with  slightly  modified  Chueh-Prausni  tz^  Henry's  and 
Fitting  programs.  The  resulting  parameters  are  given  in  Table  8. 

68 

The  low  temperature  i -butane-ethane  data  of  Skripa  et  al 
were  processed  with  the  Chueh-Prausni tz  Symfit  program.  All  of  the 
a-|2  parameters  obtained  from  this  data  were  negative  suggesting  a 
possible  thermodynamic  inconsistency. 

Prediction  of  the  Equilibrium  Data 

The  K-factor  data  were  predicted  with  the  Chueh-Prausni tz^° 
method  using  the  above  tabulated  parameters  and  also  with  the  NGPA 
version  of  the  Chao-Seader ^  method.  The  predictions  for  the  i-butane- 
carbon  dioxide  system  and  the  i -butane-ethane  system  are  given  along 
with  the  smoothed  experimental  equilibrium  constants  in  Tables  9  and 
10,  respectively. 

Thermodynamic  Consistency  Test 

The  smoothed  experimental  data  were  tested  for  thermodynamic 

consistency  using  a  method  presented  by  Chueh,  Muirbrook,  and 
72 

Prausnitz  .  The  test  is  based  on  the  Gibbs-Duhem  equation  which  for 
a  binary  system  at  constant  temperature  is 

X-,  d  In  f]  +  x2  d  In  f2  =  VL  dP/RT  10.01 

By  the  appropriate  substitutions  and  integration  the  following  point 
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by  point  consistency  test  was  arrived  at. 

Area  I  +  Area  II  +  Area  III  = 


In  <}>.|  P 


[ln  K1  +  7 77T  +  x2(ln  ^  +  ln  x-i 


10.02 


1  1  1 


where 

Area  I 


rX 

J  0 


2 


In  Kz/K^ 


10.03 


Area  II 


rX 


2 


0 


1  n  ^  /  ^1 


10.04 


Area  III 


rP  at  x? 

VL  dP 
0 


10.05 


The  three  areas  were  found  numerically  using  the  method  of 
assigned  sample  points  at  unequal  intervals  .  A  comparison  of  the 
sum  of  the  three  areas,  the  left  hand  side  of  the  equation,  with  the 
right  hand  side  of  the  equation  for  the  i-butane-carbon  dioxide  system 
and  the  i -butane-ethane  system  are  given  in  Tables  11  and  12 
respecti vely . 

The  fugacity  coefficients,  <|>..,  were  calculated  using  the 
Redlich-Kwong  equation  of  state  with  binary  interaction  parameters  and 
in  view  of  this  uncertainty  Chueh  et  al  judged  that  a  disagreement  of 
less  than  about  five  percent  constituted  very  good  thermodynamic 
consistency. 
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The  data  at  220°F  and  below  are  certainly  consistent  and  it 
is  felt  that  the  large  inconsistency  of  the  250° F  data  may  be  due  in 
part  to  the  inability  of  the  equation  of  state  to  predict  fugacities 
when  the  solvent  component  is  in  a  near-critical  situation,  and  in 
part  to  the  data  themselves. 
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Table  8.  Chueh-Prausni tz  Parameters  for  the 
i -Butane-Carbon  Dioxide  and  i -Butane-Ethane  Systems 


System 

Temperature 

°F 

hp=° 

psia 

a12 

lb.-mole/cu.ft. 

C4H10"C02 

100 

1105 

0.3267 

.  =  0.14 

J 

160 

1380 

0.5969 

220 

1450 

1.2807 

250 

1030 

0.9542 

C4H10'C2H6 

100 

463 

0.0342 

,  =  0.00 
\J 

160 

695 

0.1000 

220 

938 

0.3220 

250 

940 

0.7519 
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Table  9.  Comparison  of  the  Chueh-Prausni tz  and  Chao-Seader 


Predictions  of  the  i -Butane-Carbon  Dioxide  System 

100. 2°F 


Pressure 

i 

K 

-Butane 

%  Deviation 

Carbon 

K 

Dioxi  de 
%  Deviation 

psia 

exp 

C-S 

C-P 

exp 

C-S 

C-P 

200 

0.4479 

-4.1 

-0.3 

5.607 

-  1.6 

2.0 

300 

0.3441 

-1.8 

0.2 

3.708 

-  0.2 

5.1 

400 

0.3001 

0.1 

0.6 

2.739 

2.3 

8.2 

500 

0.2896 

-0.5 

-1.1 

2.149 

5.5 

11.3 

700 

0.3333 

-3.7 

-4.5 

1.471 

13.3 

16.2 

900 

0.4817 

- 

- 

1.145 

- 

- 

1000 

0.6164 

— 

— 

1.066 

_ 

Pressure 

i 

K  o 

160. 0°F 

-Butane 

%  Deviation 

Carbon 

K 

Dioxide 
%  Deviation 

psia 

exp 

C-S 

C-P 

exp 

C-S 

C-P 

200 

0.8659 

-4.4 

-1.2 

6.696 

10.1 

3.4 

300 

0.6579 

-3.2 

-0.7 

4.545 

8.0 

4.2 

400 

0.5523 

-0.4 

1 .2 

3.557 

3.0 

1.5 

500 

0.5083 

0.4 

1.0 

2.817 

3.5 

3.1 

700 

0.4864 

- 

4.0 

2.006 

- 

2.0 

900 

0.5341 

- 

- 

1.521 

- 

- 

1000 

0.6378 

- 

- 

1.283 

— 

— 
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Table  9.  (continued) 
220. 0°F 


i  - 

■Butane 

Carbon 

Dioxi  de 

Pressure 

psia 

Kexp 

%  Deviation 

C-S  C-P 

Kexp 

%  Deviation 
C-S  C-P 

400 

0.8642 

0.5 

1.3 

3.580 

-  5.3 

0.9 

500 

0.7752 

1.6 

3.1 

2.896 

-  7.3 

-  1.1 

600 

0.7408 

0.5 

3.3 

2.371 

-  6.6 

-  2.1 

700 

0.7497 

-3.3 

2.4 

1 .925 

-  2.2 

-  2.2 

800 

0.7978 

-9.8 

- 

1.552 

5.9 

- 

900 

0.9115 

- 

- 

1.168 

- 

- 

250. 0°F 

i -Butane 

Carbon 

Dioxide 

Pressure 

psia 

K  %  Deviation 

exp  C-S  C-P 

Kexp 

%  Deviation 
C-S  C-P 

450 

0.9756  -0.4 

0.5 

2.600 

-  5.3  -  8.3 

500 

0.9278 

1.6 

2.568 

-12.3 

600 

0.8867 

4.1 

2.020 

-10.1 

700 

0.9062 

- 

1.500 

—  - 
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Table  10.  Comparison  of  the  Chueh-Prausnitz  and  Chao-Seader 
Predictions  for  the  i -Butane-Ethane  System 

100. 6°F 


Pressure 

i 

K  „ 

-Butane 

%  Deviation 

K 

Ethane 

%  Deviation 

psia 

exp 

C-S 

C-P 

exp 

C-S 

C-P 

200 

0.4384 

-  0.1 

5.8 

2.590 

2.7 

0.1 

300 

0.3429 

2.3 

9.1 

1.826 

2.9 

0.2 

400 

0.3375 

-  5.3 

3.0 

1.447 

3.4 

0.5 

500 

0.3496 

-  9.3 

2.0 

1.236 

3.2 

0.1 

600 

0.4000 

-15.6 

0.0 

1.106 

2.8 

0.2 

700 

0.5088 

-21.3 

- 

1.030 

1 .9 

- 

Pressure 

i 

K 

160. 4°F 

-Butane 

l  Deviation 

K 

Ethane 

%  Deviation 

psia 

exp 

C-S 

C-P 

exp 

C-S 

C-P 

200 

0.8337 

-  0.2 

3.2 

3.804 

2.9 

-  1.1 

300 

0.6473 

-  0.5 

3.6 

2.513 

7.9 

3.6 

400 

0.5737 

-  2.0 

3.7 

1.958 

7.5 

2.8 

500 

0.5397 

-  2.2 

6.8 

1 .636 

6.4 

0.7 

600 

0.5660 

-  7.8 

9.5 

1.385 

8.3 

-  0.5 

700 

0.6458 

-16.2 

- 

1.205 

10.3 

- 

800 

0.8473 

-31.2 

_ 

1.054 

14.3 

— 

■ 
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Table  10.  (continued) 
219. 7°F 


assure 

i 

K  vn 

-Butane 

%  Deviation 

K 

Ethane 

%  Deviation 

psia 

exp 

C-S 

C-P 

exp 

C-S 

C-P 

400 

0.8834 

-  1.5 

-0.2 

2.165 

21.9 

21.0 

500 

0.8099 

-  1.6 

4.2 

1.785 

18.9 

9.0 

600 

0.7986 

-  3.7 

6.8 

1.493 

17.9 

3.2 

700 

0.8697 

-  9.2 

1.207 

21.1 

_ 

249. 6°F 

i-Butane  Ethane 


Pressure 

psia 

Kexp 

%  Deviation 

C-S  C-P 

Kexp 

%  Deviation 
C-S  C-P 

450 

0.9757 

-  0.5 

0.6 

2.846 

-13.3 

-16.0 

500 

0.9316 

0.2 

- 

2.000 

9.4 

- 

600 

0.9366 

- 

- 

1.323 

- 

- 
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Table  11.  Thermodynamic  Consistency  of  the 
i -Butane-Carbon  Dioxide  System 


Pressure 

psia 

xco2 

mole  frac. 

100. 2°F 

R.H.S. 

L.H.S. 

%  Deviation 

300 

0.195 

0.627 

0.598 

4.8 

500 

0.382 

1  .151 

1.131 

1.7 

900 

0.782 

2.082 

2.088 

-  0.3 

1075 

0.915 

2.335 

2.319 

0.7 

160. 0°F 

Pressure 

xco„ 

R.H.S. 

L.H.S. 

%0  Deviation 

psia 

mole  frac. 

400 

0.149 

0.407 

0.405 

0.5 

600 

0.275 

0.732 

0.734 

-  0.4 

800 

0.403 

1.024 

1.051 

-  2.6 

1082 

0.678 

1.556 

1.608 

-  3.2 

220. 0°F 

Pressure 

xco? 

R.H.S. 

L.H.S. 

%  Deviation 

psia 

mole  frac. 

600 

0.159 

0.347 

0.381 

-  9.3 

800 

0.268 

0.588 

0.626 

-  6.2 

918 

0.380 

0.793 

0.832 

-  4.8 

250. 0°F 

Pressure 

xco2 

R.H.S. 

L.H.S. 

7o  Deviation 

psia 

mole  frac. 

600 

0.072 

0.184 

0.228 

-21.2 

700 

0.128 

0.302 

0.352 

-15.4 

738 

0.202 

0.393 

0.427 

-  8.3 
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Table  12.  Thermodynamic  Consistency  Test  of  the 


i -Butane-Ethane  System 


100. 6°F 

Pressure 

psia 

^2^*6 

mole  frac. 

R.H.S. 

L.H.S. 

%  Deviation 

200 

0.261 

0.540 

0.539 

0.2 

400 

0.597 

1.253 

1.248 

0.4 

600 

0.850 

1 .769 

1.760 

0.5 

740 

0.980 

2.011 

2.018 

-  0.3 

160. 4°F 

Pressure 

psia 

Ve 

mole  frac. 

R.H.S. 

L.H.S. 

%  Deviation 

400 

0.308 

0.553 

0.578 

-  4.3 

600 

0.530 

0.943 

0.980 

-  3.8 

800 

0.738 

1.281 

1.334 

-  4.1 

219. 7°F 

Pressure 

psia 

Xc2H6 
mole  frac. 

R.H.S. 

L.H.S. 

%  Deviation 

400 

0.091 

0.135 

0.140 

-  3.9 

550 

0.242 

0.353 

0.376 

-  6.4 

700 

0.386 

0.561 

0.599 

-  6.6 

249. 3°F 

Pressure 

psia 

Xc2H6 

mole  frac. 

R.H.S. 

L.H.S. 

%  Deviation 

550 

0.115 

0.136 

0.196 

-36.5 

625 

0.200 

0.249 

0.321 

-25.3 

XI  -  CONCLUSION 


This  work  has  shown  that  in  high  pressure  vapor-liquid 
equilibrium  studies  the  densities  of  the  coexisting  phases  can  be 
obtained  conveniently  and  with  reasonable  accuracy  by  measuring  their 
refractive  indices  and  compositions.  The  designed  equilibrium  cell 
proved  to  be  convenient  to  use  and  easy  to  assemble.  The  measure¬ 
ment  of  the  refractive  indices  was  simple  and  rapid  and  the  measure¬ 
ment  of  the  compositions  was  routine. 

Although  the  refractive  indices  of  the  phases  can  be 
measured  with  high  accuracy,  there  are  some  limitations  to  the  accuracy 
of  the  obtained  densities.  These  limitations  are  caused  by  the  varia¬ 
tion  of  the  Lorentz-Lorenz  molar  refractivity  with  density  for  both 
pure  components  and  mixtures  and  also  by  the  existence  of  non¬ 
additive  terms  which  are  ignored  when  calculating  mixture  refract- 
ivities.  In  general,  these  limitations  seldom  cause  inaccuracies 
greater  than  one  or  two  percent  in  the  calculated  density.  The  experi¬ 
mental  uncertainty  in  the  composition  of  the  phases  can  also  lead  to 
slight  inaccuracies  in  the  calculated  refracti vities . 

When  data  of  higher  precision  are  required  the  nature  and 
magnitude  of  the  variation  of  the  refractivity  with  density  and  com¬ 
position  must  be  considered.  This  information  can  best  be  obtained 
by  measuring  the  refractive  index  and  density  in  single  phase  regions 
where  it  is  much  less  difficult  to  obtain  accurate  volumetric  data. 

The  information  obtained  will  prove  equally  valid  in  the  two  phase 
region  and  improve  the  accuracy  of  the  obtained  densities. 

The  use  of  refractive  index,  therefore,  has  the  possibility 
of  becoming  a  very  powerful  aid  in  the  experimental  determination  of 
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densities  in  high  pressure  two  phase  systems  and  in  single  phase  systems 
in  the  near  critical  regions. 


NOMENCLATURE 


d  distance  between  two  plates  of  an  interferometer 

f  fugacity 

k  number  of  interference  fringe  shifts 

i  angle  of  incidence 

m  induced  electrical  moment 

n  refractive  index 

r  angle  of  refraction 

x  composition,  mole  fraction 

a  polarizability 

van  Laar  interaction  parameter 
3  angle  of  refraction  in  the  window 

y  minimum  deviation  angle  measured  by  autocol limation 

6  prism  angle 

e  permittivity 

p  density 

A  wavelength 

<t>  fugacity  coefficient 

A  prism  angle 

B  critical  angle  of  refraction 

B  second  refractivity  virial  coefficient 

C  third  refractivity  virial  coefficient 

D  minimum  deviation  angle 

F  electric  field  strength 

H  Henry's  law  constant 

K  vapor-liquid  equilibrium  constant 
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M  molecular  weight 

N  Avogadro's  number 

P  pressure 

Lorentz-Lorenz  refractivity 

zero  density  limit  of  the  Lorentz-Lorenz  molar  refractivity 
T  temperature 

V  molar  volume 
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APPENDIX  A 


DERIVATION  OF  THE  LORENTZ-LORENZ  EQUATION 

In  a  parallel  plate  condenser  across  which  a  potential 
difference,  V,  is  maintained,  the  field,  F,  acting  on  a  single  mole¬ 
cule  in  the  dielectric  can  be  calculated  by  placing  an  imaginary 
sphere  around  the  molecule.  The  radius  being  small  compared  with 
the  dimensions  of  the  condenser  but  large  compared  with  the  molecular 
dimensions.  Outside  the  sphere  the  material  is  treated  as  homogeneous 
while  inside  the  sphere  it  is  treated  as  if  it  were  made  up  of 
individual  molecules. 

The  field  acting  on  the  molecule  is  due  to  three  sources. 

1.  The  charge  on  the  condenser  plates  and  the  dielectric  surfaces 
adjacent  to  them  giving  rise  to  a  field  equal  to  the  applied 
field 


F,  -  M  .  E 
I  £ 


A-l 


2.  The  polarization  charges  on  the  spherical  surface  around  the 
molecule 


n  Surface  Charqe  Density  •  cos  0  •  d  Area 

2 

0  e(Distance) 


-TT  ? 

r  P  cos  0  •  cos  0  •  2nr  sin  0  d  0 


e  r 


F2  - 


n 


P  cos  0  sin©  d© 


A-2 
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A-3 


3.  The  molecules  inside  the  spherical  region  assumed  to  be  in  a 

spherically  symmetric  arrangement  thereby  providing  no  net  field. 


is 


A-4 


Therefore,  the  field  acting  on  a  molecule  in  the  dielectric 


F  =  F1  +  F2  +  F3  =  E  +  TT  A'5 

The  permittivity  or  dielectric  constant,  eQ,  of  an  isotropic 
material  is  related  to  the  polarization,  P,  the  product  of  the  average 
induced  moment  per  molecule  produced  by  an  applied  field,  E,  and  the 
number  of  molecules  per  unit  volume,  by 

P  =  In  E(eo  ‘  1)  A'6 

where  e  is  the  permittivity  of  free  space  in  rationalized  units. 
Substituting  into  Equation  A-5  for  P  and  simplifying  we  have 

F  =  (eQ  +  2)  E/3  A-7 

By  definition,  the  molar  polarization  at  visible  frequencies 
is 


P  =  N,  a  F 
1  e 

and  the  number  of  molecules  per  unit  volume,  N-j ,  is 


A-8 


Np/M 


A-9 
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Combining  Equations  A-6,  A-7,  A-8,  and  A-9, 
and  simplifying  we  have  the  Cl ausi us-Mossi tti  equation  applicable  at 
visible  frequencies. 


£o  ~  1  M  .  4nN 
eQ  +  2  p  3  e  ae 


A- 1 0 


From  Maxwell's  electromagnetic  theory  it  follows  that  the  permittivity 

2 

of  a  material,  eQ,  is  equal  to  the  square  of  the  refractive  index,  n  , 
measured  at  the  same  frequency.  Therefore  we  have  the  Lorentz-Lorenz 
equati on 


n2  -  1  M  _  4nN 
n2  +  2  p  "  3  e 


A-l  1 


APPENDIX  B 


COMPARISON  OF  THE  CHAO-SEADER  AND  CHOEH-PRAUSNITZ 
APPROACHES  TO  VAPOR-LIQUID  EQUILIBRIUM  PREDICTIONS 

Computers  have  significantly  reduced  the  time  and  effort 
required  to  calculate  the  vapor-liquid  equilibrium  data  needed  in  the 
design  of  contacting  equipment.  In  general,  the  prediction  methods 
are  semi -empiri cal ,  being  based  on  a  thermodynamic  foundation  with  a 
significant  amount  of  experimental  data  for  the  evaluation  of  para¬ 
meters  . 

Two  such  methods,  one  recent,  the  Chueh-Prausnitz\  and 

2  3  4  5 

the  other  relatively  old,  the  Chao-Seader  *  *  *  ,  have  a  similar 
thermodynamic  foundation  but  differ  mainly  in  the  amount  of  experi¬ 
mental  data  used  in  the  evaluation  of  parameters. 

The  Thermodynamic  Method 

The  thermodynamic  method  has  been  presented  in  detail  by 

g 

Prausnitz  and  only  an  outline  of  its  basic  features  will  be  given 
here.  The  thermodynamic  criteria  which  determine  the  state  of  equili¬ 
brium  between  two  phases  of  a  multicomponent  system  are  equality  of 
temperature  and  pressure  in  both  phases  and  for  each  component 
equality  of  its  fugacity  in  both  phases.  In  order  to  find  the  condi¬ 
tions  which  satisfy  these  criteria,  it  is  necessary  to  have  some 
method  for  evaluating  the  fugacity  of  each  component  in  each  phase. 

The  fugacity  of  any  component  can  be  obtained  from  the 
volumetric  data  for  the  mixture  in  question  by  integrating  according 
to  Equation  B-l  from  a  pressure  of  zero  to  the  operating  pressure. 

This  equation  applies  to  both  the  vapor  and  the  liquid  phases.  For 
the  liquid  phase,  x  replaces  y. 
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In  f.  =  In 


RT 


[(!£) 


an. 'T,P,n. 


-  ]  dp 


B-l 


This  relatively  simple  equation  involves  a  presupposed  know¬ 
ledge  of  the  behavior  at  constant  temperature  of  the  total  volume  of 
the  mixture,  with  respect  to  a  slight  compositional  change, from  zero 
to  the  operating  pressure.  For  the  vapor  phase  this  knowledge  can  be 
provided  by  an  equation  of  state.  However,  for  the  liquid  phase  the 
necessary  density  range  over  which  the  integration  is  to  be  performed 
is  usually  beyond  the  scope  of  any  equation  of  state.  Some  success 
with  the  Redlich-Kwong  and  Benedict-Webb-Rubin  equations  have  been 
reported^*^’^. 

An  alternative  method  for  obtaining  liquid  phase  fugacities 
is  by  defining  an  ideal  solution  and  using  a  solutions  theory  to 
determine  deviations  from  this  ideal  behavior.  This  technique  has 
been  used  in  both  the  Chao-Seader  and  Chueh-Prausnitz  methods  to  deter¬ 
mine  the  liquid  phase  fugacities.  Also,  both  methods  determine  vapor 
phase  fugacities  using  a  form  of  the  Redlich-Kwong  equation  of  state. 

Once  methods  for  calculating  liquid  and  vapor  fugacities  are 
available,  an  iterative  procedure  can  be  set  up  to  satisfy  the  thermo¬ 
dynamic  criteria  with  the  restriction  that  the  liquid  and  vapor  com¬ 
positions  both  sum  to  1.000. 


Vapor  Phase  Fugacities 

In  both  the  Chao-Seader  and  the  Chueh-Prausnitz  methods  the 
vapor  phase  fugacities  are  calculated  with  the  Redlich-Kwong  equation 


of  state. 
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P  V  _  V  a 

RT  -  V-b  -  rj3/2  (v+b) 

A  vapor-phase  fugacity  coefficient,  ,  is  defined  by  rewriting 
Equation  B-l  as  follows 


In  f^v  =  In  y.j  P  +  In  <f>. 


B-3 


Performing  the  required  integration  on  the  R-K  equation  gives 


In  <j>.  = 


V  +  bi 


ai 


V-b  V-b 


RT 


_  ln  rV+bx 

2/3  k  n  V 


B-4 


To  calculate  <j>.  a  value  of  V,  the  total  volume,  is  obtained 
by  solving  the  R-K  equation  at  a  given  T,  P,and  y.  This  value  is 
then  substituted  into  equation  B-4  adhering  to  the  mixing  rules  for 
the  constants  as  given  below. 


a 


z  y 


1  1 


b 


z  y. 


bi 


2  2  5 

0.4278  R^  T/*° 
c  • 


0.0867 


B-5 


B-6 


Chueh  and  Prausnitz  modified  the  mixing  rules  to  include  a  binary 
interaction  parameter, k. . . ,  and  substituted  specific  component  para- 

*  J 

meters  for  the  universal  constants.  Their  mixing  rules  are  given  in 
equations  B-7  to  B-l 3 

(n  +  n  )  R2t  2-5 

r  r  3i  ai  °ii 

=  n  yi  yj  ■iJ  .  a1d  -  1  2  -p— J- 

1  J  ci  j 


a 


B-7 
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ft  b  # 

l  yi  bi  ’  bi  P  T 

ci  ci 

b  = 

B-8 

ai  = 

2  l  h  aki 

B-9 

where 

Sj  RTCii 

P 

Cij 

-  _ 1 J  1 J 

vc.. 

10 

B-10 

vc.. 

10 

=  1/2(V1/3  +  v1/3)3 

Ci  cj 

w.  +  w. 

=  0.291  -  0.08  (-1-  -2-  ^ 

B-l  1 

ZC.. 

10 

B-12 

TC.. 

1J 

B-l  3 

Liquid  Phase  Fugacities 

An  ideal  solution  can  be  defined  as  one  in  which  the  fugacity 
of  a  component  in  the  solution  is  equal  to  the  product  of  its  mole 
fraction  and  the  fugacity  of  the  pure  component  at  the  same  temperature 
and  pressure  over  the  entire  composition  range  and  a  non-zero  range  of 
temperature  and  pressure. 

Chao-Seader  Pure  Component  Fugacities 

The  pure  component  fugacity  in  the  Chao-Seader  correlation 
is  obtained  by  multiplying  the  operating  pressure  by  a  fugacity 
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coefficient,  v?.  This  fugacity  coefficient  is  calculated  from  a 
corresponding  states  correlation  in  Pitzer's  modified  form, 

In  v9  =  in  »  +  w.  In  v^1 )  b-14 

where  and  '  are  functions  of  the  reduced  temperature  and 
pressure  only. 


Chao-Seader  Deviation  From  Ideal  Solution  Behavior 

The  deviation  from  ideal  solution  behavior  is  represented  by 
the  activity  coefficient, y- , obtained  from  Hildebrand's^  regular 
sol ution  model . 


In 


Yi 


Vi(«i  -  <5)2 

RT 


B  - 1 5 


where  V..  is  the  pure  liquid  molar  volume  at  25°C  and  its  vapor  pressure 
and  6  is  free  energy  density  function  defined  by 


a  -  C^)1/2  =  [^]1/2 


B  - 1 6 


where  Ay  is  the  latent  heat  of  vaporization  at  25°C. 

The  total  solution  solubility  parameter,  T,  is  defined  by 


For  components  which  are  super-critical  or  near  critical  at  25°C  such 
as  methane,  ethane  and  propane  the  6^  and  V..  parameters  can  be  evaluated 
from  a  best  fit  of  binary  vapor-liquid  equilibrium  data  for  the  light 
component  with  several  of  the  heavier  components. 

All  the  necessary  terms  are  defined  and  the  liquid  phase 
fugacity  is  given  by  equation  B-18. 
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In  f\  =  In  (x.  P)  +  In  v?  +  In  Y° 


B-18 


Chueh-Prausni tz  Pure  Component  Fugacities 

In  the  Chueh-Prausni tz  method  a  distinction  is  made  between 
solvents  and  solutes,  a  solute  being  defined  as  any  component  whose 
reduced  temperature  is  greater  than  0.93,  and  a  solvent  being  any 
component  whose  reduced  temperature  is  less  than  0.93. 

A  solvent  pure  component  fugacity  is  obtained  by  calculating 
the  fugacity  of  the  vapor  at  the  solvent's  vapor  pressure  with  the  R-K 
equation  of  state  and  reducing  it  to  zero  pressure  by  a  Poynting 
correction  according  to  equation  B-19. 


In  f°,P=0 


In  f?’P=P: 


B-19 


This  correction  is  necessarily  hypothetical  because  the 
saturated  liquid  molar  volume,  V^,  is  assumed  constant  from  the  satura¬ 
tion  pressure,  P  ,  to  zero  pressure  and  in  fact  no  liquid  can  exist 
below  its  saturation  pressure. 

A  solute  pure  component  fugacity  is  calculated  from  its 
Henry's  Law  constant  in  a  solvent,  and  is  also  reduced  to  zero  pressure 
by  a  Poynting  correction 

V  p=ps  •  p 

o,P=Q  o,P=Ps.  VL.  SJ  Ks. 

in  f1d  -  in  Htj  3  (-Hit — i)  B-20 

where  V,  is  the  partial  molar  volume  of  the  solute  i  in  the  solvent  j 
Li 

at  infinite  dilution.  The  use  of  Henry's  Law  constants  has  two  advan¬ 
tages  over  'hypothetical  liquid'  fugacities.  It  is  an  experimentally 
measurable  quantity  and  the  pure  solute  fugacity  is  a  function  of  the 
solvent  in  which  it  is  dissolved.  However,  the  Poynting  correction  is 
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still  hypothetical. 

Chueh-Prausni tz  Deviations  From  Ideal  Solution  Behavior 

The  deviations  from  ideal  solution  behavior  are  calculated 
using  the  van  Laar  equation  with  a  dilation  parameter,  n. 

In  y-  °  =  a^.  Vc  q?  [1  +  3n ^  q^]  B-21 

where 


n-  =  dilation  parameter  for  solutes  in  the  form  of  an  empirical 
function  of  T  but  is  set  equal  to  zero  for  solvents. 


the  van  Laar  interaction  constant  which  must  be  obtained 
from  binary  vapor-liquid  equilibrium  data. 


Partial  Molar  Volumes 

The  partial  molar  volumes  in  all  of  the  Poynting  corrections 
are  calculated  using  the  modified  Redl ich-Kwong  equation  of  state. 


RT 

V-b 


(1  + 


RT 

(V-b)2 


V-b 


)  - 


^n 


n  ab. 

2U1  xiaki}  '  TvTb) 

V( V+b)  T1/2 

2V+b 

V2( V+b)2 


B-22 


The  liquid  molar  volume,  V,  is  not,  however,  obtained  from  the  equation 

of  state  but  rather  from  the  corresponding  states  correlation  of 

1 2 

Lyckman  and  Eckert  .  The  mixing  rules  are  the  same  as  those  used  for 


the  vapor  phase  fugacity. 
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Having  defined  all  of  the  necessary  terms  the  liquid  phase 
fugacities  can  be  written  according  to  equation  B -23 


In  f. 


i  t  *  05P=0,  .  ,  P=0  A  ViP 

ln(x.f.  )  +  In  y1  + 


B-23 


Chueh-Prausni tz  Program 

For  convenient  use  in  the  program  the  following  properties 
are  fitted  to  polynomials  in  temperature: 

1.  The  pure  component  zero  pressure  solvent  fugacities. 

2.  The  Henry's  law  constants  reduced  to  zero  pressure  as 
determined  from  solubility  data  are  fitted  to  a  poly¬ 
nomial  for  each  binary  solvent-solute  pair. 

3.  The  zero  pressure  van  Laar  interaction  constants  as 
determined  from  binary  vapor-liquid  equilibrium  data 
are  fitted  to  a  polynomial  for  each  solvent-solvent 
pair  and  solvent-solute  pair. 


Nomenclature 


f 

y 

p 

R 

T 

V 


n 

a  j  b  j 


Y 

6 


H 


a 

q 

n 

IX) 


fugaci ty 

vapor  phase  mole  fraction 

pressure 

gas  constant 

absolute  temperature 

molar  volume 

number  of  moles 

Redl i ch-Kwong  constants 

vapor  phase  fugaci ty  coefficient 

critical  temperature 

critical  pressure 

critical  volume 

critical  compressibility  factor 

binary  interaction  parameter 

pure  component  fugaci ty  coefficient 

activity  coefficient 

solubility  parameter 

heat  of  vaporization 

free  energy  of  vaporization 

Henry's  law  constant 

van  Laar  interaction  parameter 

critical  volume  fraction 

dilation  parameter 

accentric  factor 
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APPENDIX  C 


EQUILIBRIUM  CELL  DIMENSIONS 
NOTE:  ALL  DIMENSIONS  ARE  IN  INCHES 


90 


91 


FIG.  24  DIMENSIONS  OF  THE  CELL  MAIN  BODY 
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FIG.  25  CROSS-SECTIONS  OF  THE  CELL  MAIN  BODY 
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FIG.  26  DIMENSIONS  OF  THE  END  CYLINDERS 
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FIG.  27  DIMENSIONS  OF  THE  PISTONS 


APPENDIX  D 


EXPERIMENTAL  DATA 


This  Appendix  lists  the  experimental  data  points  taken  during 
the  course  of  this  work.  The  actual  experimental  data  on  the  pure  sub¬ 
stances  consists  of  three  parameters  only:  the  temperature,  the  pres¬ 
sure,  and  the  refractive  index  at  the  wavelength  6328A- 


In  Table  13:  the  figures  in  the  column  headed  by  'Specific 
Volume  Exp.1  are  interpolated  data  taken  from  Sage,  Webster  and  Lacey 
(Ref.  67  main  text); 

the  figures  in  the  column  headed  by  'Molar 
Refractivity  Exp.'  are  the  values  of  the  Lorentz-Lorenz  molar  refracti vi ty , 


R 


LL 


n2-l  m 
n2+2  p 


the  refractive  index  is  the  experimental  data  point  and  the  density  is 
that  obtained  from  Sage,  Webster  and  Lacey; 

the  figures  in  the  column  headed  by  'Molar 
Refractivity  Calc.'  are  the  values  obtained  from  the  virial  refractivity 
expansion  where 

Rll  =  Rll  +  B/V  +  C/V2 
o 


V  has  units  of  cc/g-mole.  The  values  of  Band  C  were  calculated  so  as 
to  give  a  best  fit  in  the  least  squares  sense  of  the  data  in  the  column 
headed  'Molar  Refractivity  Exp.'  and  are  given  in  Appendix  E; 

the  figures  in  the  column  headed  'Specific  Volume 
Calc.'  were  obtained  by  substituting  the  experimental  refractive  index 
into  the  values  listed  in  the  column  'Molar  Refractivity  Calc.'. 


In  Table  14:  the  headings  have  a  similar  connotation  to  those 
in  Table  13  except  that  the  figures  in  the  column  headed  by  'Specific 
Volume  Exp.'  were  interpolated  from  the  volumetric  data  of  Sage  and 
Lacey  (Reference  69  main  text). 

In  Tables  15  and  16  the  figures  in  the  column  headed  by  'RLL' 
were  obtained  from  the  virial  expansion  using  the  coefficients  listed  in 
Appendix  E  and  the  figures  in  the  column  headed  by  'Spec.  Vol . '  have  the 
same  connotation  as  those  in  the  column  headed  'Specific  Volume  Calc.' 
in  Tables  14  and  15. 


In  Tables  17  and  18  the  actual  experimental  data  consists  of 
the  temperature,  the  pressure,  and  the  liquid  and  vapor  compositions  and 
refractive  indices.  The  figures  listed  under  the  column  headed  by  V* 
were  calculated  as  outlined  on  page  37  of  the  main  text. 
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TABLE  14.  EXPERIMENTAL  REFRACTIVE  INDEX  DATA  FOR  I  “-BUTANE  LIQUID 
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TABLE  15.  EXPERIMENTAL  REFRACTIVE  INDEX  DATA  FOR  CARBON  DIOXIDE 


TEMPERATURE 
DEG.  F 

PRESSURE 

PS  I  A 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU . FT. /LB 

100.17 

14.9 

1.00040 

6.645 

9.85785 

100.17 

97.0 

1.00263 

6.645 

1.51382 

100.17 

168.6 

1.00468 

6.645 

0.85285 

100.17 

244.7 

1.00697 

6.645 

0.57246 

100.17 

321.2 

1.00942 

6.645 

0.42408 

100.13 

395.4 

1.01197 

6.646 

0.33408 

100.17 

474.2 

1.01478 

6.646 

0.27068 

99.96 

547.7 

1.01764 

6.646 

0.22682 

99.96 

694.3 

1.02427 

6.646 

0.16508 

100.10 

767.0 

1.02806 

6.646 

0.14291 

100.17 

818.9 

1.03104 

6.646 

0.12926 

100.17 

883.5 

1.03518 

6.646 

0.11412 

100.17 

933.0 

1.03894 

6.645 

0.10318 

100.17 

1011.1 

1.04694 

6 .645 

0.08572 

100.17 

1078.8 

1.05364 

6.644 

0.07509 

100.00 

1133.8 

1.06275 

6.642 

0.06428 

100.00 

1177.8 

1.07331 

6.640 

0.05512 

100.00 

1202.1 

1.08319 

6.637 

0.04865 

100.06 

1224.5 

1.09438 

6.633 

0.04296 

100.10 

1245.4 

1.10738 

6.628 

0.03783 
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TABLE  15.  (CONTINUED) 

TEMPERATURE 
DEG.  F 

PRESSURE 

PSIA 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU.FT./LB 

100.17 

1263.6 

1.11778 

6.623 

0.03454 

100.03 

1293.1 

1.13061 

6.617 

0.03120 

99.86 

1308.8 

1.13541 

6.614 

0.03012 

100.00 

1343.5 

1.14212 

6.610 

0.02872 

99.93 

1370.3 

1.14685 

6.607 

0.02781 

100.10 

1403.2 

1.15061 

6.605 

0.02713 

100.10 

1443.5 

1.15431 

6.603 

0.02649 

100.20 

1470.3 

1.15654 

6.601 

0.02612 

100.27 

1476.7 

1.15667 

6.601 

0.02610 

100.17 

1480.4 

1.15706 

6.601 

0.02604 

TABLE  15.  (CONTINUED) 


TEMPERATURE 
DEG.  F 

PRESSURE 

PS  I  A 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU.FT./LB 

100.06 

139.5 

1.00382 

6.645 

1.04463 

100.00 

262.2 

1.00737 

6.645 

0.54211 

100.00 

429.0 

1.01294 

6.646 

0.30896 

100.17 

627.5 

1.02081 

6.646 

0.19248 

99.86 

719.7 

1.02540 

6 .646 

0.15777 

99.75 

812.2 

1.03041 

6 .646 

0.13193 

99.79 

940.9 

1.03932 

6.645 

0.10218 

99.82 

1093.1 

1.05569 

6.643 

0.07235 

99.82 

1162.9 

1.06938 

6.641 

0.05821 

99.93 

1192.8 

1.07883 

6.638 

0.05130 

100.03 

1214.0 

1.08816 

6.635 

0.04594 

c— • 

• 

o 
o 
«■— • 

1236.8 

1.10160 

6.630 

0.03995 

100.41 

1267.0 

1.11757 

6.623 

0.03460 

100.41 

1299.5 

1.13059 

6.617 

0.03121 

100.03 

1359.9 

1.14463 

6.609 

0.02823 

99.82 

1350.2 

1.14426 

6.609 

0.02830 

99.65 

1412.9 

1.15241 

6.604 

0.02682 

100.17 

1488.6 

1.15797 

6.600 

0.02589 

100.34 

1499.1 

1.15817 

6.600 

0.02586 
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TABLE  15.  (CONTINUED) 


TEMPERATURE 
DEG.  F 

PRESSURE 

PS  I  A 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL. 
CU.FT./LB. 

159.89 

63.1 

1.00160 

6.645 

2.48115 

159.89 

178.9 

1.00450 

6.645 

0.88714 

159.93 

294.4 

1.00749 

6.645 

0.53334 

159.93 

411.4 

1.01074 

6.646 

0.37205 

160.03 

632.6 

1.01759 

6.646 

0.22754 

159.93 

743.7 

1.02150 

6.646 

0.18631 

160.03 

839.8 

1.02503 

6.646 

0.16009 

160.03 

949.7 

1.02938 

6 .646 

0.13652 

160.06 

1051.1 

1.03372 

6.646 

0.11905 

160.10 

1175.9 

1.03960 

6.645 

0.10146 

160.13 

1294.3 

1.04582 

6.645 

0.08779 

160.10 

1420.5 

1.05331 

6.644 

0.07556 

160.10 

1421.3 

1.05336 

6.644 

0.07549 

160.06 

1480.3 

1.05721 

6.643 

0.07044 

160.06 

1482.0 

1.05748 

6.643 

0.07012 

TABLE  15.  (CONTINUED) 


TEMPERATURE 
DEG.  F 

PRESSURE 

PS  I  A 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU. FT. /LB 

219.72 

89.4 

1.00200 

6.645 

1.98643 

219.72 

126.3 

1.00286 

6.645 

1.39487 

219.68 

163.3 

1.00369 

6.645 

1.07980 

219.68 

200.2 

1.00456 

6.645 

0.87449 

219.72 

237.2 

1.00541 

6.645 

0.73712 

219.72 

274.1 

1.00627 

6.645 

0.63710 

219.65 

311.1 

1.00715 

6.645 

0.55840 

219.72 

348.0 

1.00804 

6.645 

0.49705 

219.72 

385.0 

1.00900 

6.645 

0.44372 

219.72 

421.9 

1.00985 

6.645 

0.40550 

219.72 

458.9 

1.01079 

6.645 

0.37047 

219.65 

495.8 

1.01172 

6.646 

0.34104 

219.65 

532.8 

1.01267 

6 . 646 

0.31554 

219.65 

569.7 

1.01369 

6.646 

0.29217 

219.65 

606.7 

1.01467 

6 . 646 

0.27266 

219.65 

643.6 

1.01562 

6 . 646 

0.25614 

219.62 

680.6 

1.01662 

6 . 646 

0.24080 

219.72 

717.5 

1.01762 

6.646 

0.22720 

219.72 

754.5 

1.01870 

6.646 

0.21410 

219.72 

791.4 

1.01978 

6.646 

0.20245 

219.72 

828.4 

1.02081 

6.646 

0.19246 

219.68 

852.4 

1.02159 

6.646 

0.18553 
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TABLE  15.  (CONTINUED) 
TEMPERATURE  PRESSURE 


N 


RLL 


SPEC.  VOL. 


DEG.  F 

PSIA 

CC/G-MOLE 

CU.FT./LB 

219.68 

865.3 

1.02189 

6.650 

0.18312 

219.68 

902.3 

1.02313 

6.650 

0.17331 

219.79 

976.2 

1.02526 

6.651 

0.15880 

219.79 

939.2 

1.02429 

6.651 

0.16506 

219.82 

1013.1 

1.02643 

6.651 

0.15177 

219.82 

1050.1 

1.02764 

6.651 

0.14516 

219.82 

1087.0 

1.02872 

6.651 

0.13975 

219.93 

1124.0 

1.02998 

6.651 

0.13391 

219.86 

1160.9 

1.03112 

6.652 

0.12903 

219.89 

1197.9 

1.03241 

6.652 

0.12392 

219.86 

1234.8 

1.03364 

6.652 

0.11944 

219.93 

1271.8 

1.03494 

6.652 

0.11500 

219.89 

1308.7 

1.03617 

6.652 

0.11114 

219.93 

1327.2 

1.03688 

6.652 

0.10902 

219.96 

1345.7 

1.03749 

6.652 

0.10725 

219.96 

1364.1 

1.03817 

6.652 

0.10536 

219.96 

1382.6 

1.03886 

6.652 

0.10350 

219.96 

1401.1 

1.03950 

6.652 

0.10182 

220.00 

1419.6 

1.04021 

6.652 

0.10004 

220.00 

1438.0 

1.04092 

6.653 

0.09832 

220.00 

1456.5 

1.04150 

6.653 

0.09696 

220.00 

1475.0 

1.04221 

6.653 

0.09535 
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TABLE  15.  (CONTINUED) 


TEMPERATURE 
DEG.  F 

PRESSURE 

PS  I  A 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU.FT./LB 

249.82 

232.9 

1.00500 

6.645 

0.79869 

249.79 

340.5 

1.00743 

6.645 

0.53719 

249.79 

456.6 

le01007 

6.645 

0.39683 

249.86 

633.9 

1.01433 

6.646 

0.27900 

249.82 

895.3 

1.02112 

6.646 

0.18962 

249.82 

761.1 

1.01756 

6.646 

0.22786 

249.86 

868.0 

1.02046 

6.646 

0.19575 

249.93 

980.4 

1.02349 

6 . 646 

0.17053 

249.96 

1087.0 

1.02645 

6.646 

0.15157 

249.96 

1183.9 

1.02931 

6.646 

0.13682 

250.20 

1337.2 

1.03396 

6.646 

0.11818 

250.20 

1444.9 

1.03732 

6.645 

0.10762 

250.20 

1580.9 

1.04173 

6.645 

0.09633 
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TABLE  16#  EXPERIMENTAL  REFRACTIVE  INDEX  DATA  FOR  ETHANE 


TEMPERATURE 
DEG  •  F 

PRESSURE 

PSIA 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU.FT./LB 

100.55 

29.7 

1.00143 

11.231 

6.26607 

100.55 

59.9 

1.00282 

11.233 

3.17966 

100.55 

152.2 

1.00746 

11.239 

1.20653 

100.55 

244.4 

1.01274 

11.246 

0.70762 

100.58 

304.6 

1.01635 

11.250 

0.55206 

100.58 

392.0 

1.02256 

11.257 

0.40082 

100.58 

487.6 

1.03051 

11.266 

0.29704 

100.58 

576.5 

1.04008 

11.275 

0.22667 

100.58 

668.8 

1.05441 

11.286 

0.16762 

100.58 

705.7 

1.06289 

11.292 

0.14533 

100.58 

733.4 

1.07192 

11.297 

0.12737 

100.58 

751.8 

1.08072 

11.302 

0.11372 

100.55 

770.3 

1.09371 

11.306 

0.09826 

100.55 

781.0 

1.10480 

11.309 

0.08808 

no 


TABLE  16.  (CONTINUED) 


TEMPERATURE 
DEG.  F 

PRESSURE 

PS  I  A 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU.FT./LB 

100.55 

788.7 

1.11696 

11.310 

0.07914 

100.27 

797.9 

1.13329 

11.309 

0.06968 

100.51 

816.4 

1 ® 14620 

11.307 

0.06370 

o 

o 

• 

853.3 

1.16277 

11.301 

0.05740 

100.41 

908.6 

1.17192 

11.296 

0 .05444 

100.41 

945.5 

1.17778 

11.293 

0.05271 

100.41 

1037.8 

1.18669 

11.288 

0.05027 

99.48 

1037.8 

1.18669 

11.288 

0.05027 

99.48 

1111.6 

1.19221 

11.284 

0.04888 

99.65 

1168.8 

1.19548 

11.281 

0.04809 

99.48 

1185.4 

1.19662 

11.280 

0.04782 

99.65 

1277.6 

1.20128 

11.277 

0.04675 

99.65 

1462.1 

1.20869 

11.270 

0.04514 

99.65 

1554.4 

1.21188 

11.268 

0.04449 

Ill 


TABLE  16.  (CONTINUED) 
TEMPERATURE  PRESSURE 


N 


RLL 


SPEC*  VOL. 


DEG.  F 

PS  I  A 

CC/G-MOLE 

CU.FT./LB 

160*44 

13*2 

1.00053 

11.230 

16.88352 

160.44 

38.3 

1.00164 

11.232 

5.45853 

160*44 

69.2 

1.00295 

11.234 

3.04059 

160*44 

84*7 

1.00356 

11.234 

2.52396 

160*44 

143.0 

1.00623 

11.238 

1,44441 

160*44 

253*5 

1.01144 

11.244 

0.78755 

160*44 

345*7 

1*01623 

11.250 

0.55626 

160.44 

437.8 

1.02142 

11.256 

0.42196 

160*44 

561.3 

1.02914 

11.264 

0.31085 

160*44 

677.4 

1.03776 

11.273 

0.24045 

160*44 

762.2 

1.04503 

11.279 

0.20202 

160*44 

854.4 

1.05412 

11.286 

0.16851 

160*44 

898.6 

1.05884 

11.290 

0.15517 

160*44 

972.3 

1.06762 

11.295 

0.13532 

160*44 

1046.0 

1.07779 

11.300 

0.11792 

160.44 

1119.8 

1.08941 

11.305 

0.10288 

160.44 

1156.6 

1.09541 

11.307 

0.09655 

160.44 

1248*8 

1.11037 

11.310 

0.08374 

160.44 

1322.5 

1.12179 

11.310 

0.07608 

160.44 

1377.8 

1.12960 

11.310 

0.07162 

160*44 

1411.7 

1.13424 

11.309 

0.06921 

160.44 

1451.5 

1.13837 

11.308 

0.06720 

160*44 

1469.9 

1.14048 

11.308 

0.06622 

160*44 

1543.6 

1.14804 

11.306 

0.06293 

TABLE  16.  (CONTINUED) 


TEMPERATURE 
DEG.  F 

PRESSURE 

PSI  A 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU.FT./LB 

219.48 

41.2 

1.00153 

11.232 

5.86831 

219.51 

115.0 

1.00439 

11.235 

2.04578 

219.65 

188.7 

1.00740 

11.239 

1.21553 

219.93 

299.3 

1.01213 

11.245 

0.74337 

219.86 

409.8 

1.01704 

11.251 

0.52992 

219.31 

520.4 

1.02254 

11.257 

0.40117 

219.44 

612.6 

1.02720 

11.262 

0.33283 

219.51 

723.1 

1.03340 

11.269 

0.27155 

219.51 

815.3 

1.03892 

11.274 

0.23339 

219.48 

903.8 

1.04452 

11.279 

0.20431 

219.48 

977.5 

1.04966 

11.283 

0.18341 

219.31 

1051.2 

1.05516 

11.287 

0.16535 

219.44 

1124.9 

1.06041 

11.291 

0.15119 

219.51 

1198.6 

1.06618 

11.294 

0.13821 

219.65 

1272.4 

1.07241 

11.298 

0.12651 

219.89 

1346.1 

1.07844 

11.301 

0.11697 

220.00 

1419.8 

1.08494 

11.303 

0.10818 

220.00 

1493.5 

1.09146 

11.306 

0.10062 

219.65 

1567.2 

1.09801 

11.307 

0.09404 
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TABLE  16.  (CONTINUED) 


TEMPERATURE  PRESSURE 


DEG.  F 

PS  I A 

249.79 

o 

• 

qo 

1 . 

249.44 

301.8 

1. 

249.55 

506.0 

1. 

249.55 

696.6 

1. 

249.55 

872.1 

1. 

249.55 

1024.7 

1. 

249.55 

1163.6 

1. 

249.55 

1293.4 

1. 

249.58 

1426.1 

1. 

249.48 

1507.5 

1. 

249.48 

1565.8 

1. 

N 

RLL 

CC/G-MOLE 

SPEC.  VOL 
CU.FT./LB 

00303 

11.234 

2.95480 

01157 

11.245 

0.77767 

02021 

11.255 

0.44636 

02933 

11.265 

0.30839 

03877 

11.274 

0.23389 

04759 

11.281 

0.19101 

05629 

11.288 

0.16184 

06493 

11.293 

0.14061 

07409 

11.298 

0.12351 

07960 

11.301 

0.11511 

08370 

11.303 

0.10957 

Table  17.  Experimental  Equilibrium  Phase  Properties  for  the  i -Butane-Carbon  Dioxide  System 
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Table  18.  Experimental  Equilibrium  Phase  Properties  for  the  i -Butane-Ethane  System 
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APPENDIX  E 


PURE  COMPONENT  REFRACTI VITI ES 


Table  19 

Pure  Component  Refractivi 

ties 

Compound 

rLL 

Temperature 

B 

C/B 

LLo 

Range 

cc/g 

-mole 

deg.  F 

cc/g-mole 

cc/g-mole 

NaD 

He-Ne 

n-Butane 

20. 621 

20.56 

100.24 

254±30 

-109 

30  -  200 

1 41  ±  4 

-  991 

i -Butane 

20. 681 

20.62 

100.0 

260±25 

-107 

159.9 

1 74±25 

-111 

220.0 

1 25±25 

-108 

250.0 

85±25 

-106 

100  -  250 

66±  6 

-109 

30  -  200 

1 43±  4 

-  691 

Propane 

15.921 

15.87 

30  -  200 

83. 5±  3 

-  771 

Ethane 

11. 261 

11.23 

30  -  200 

27.8+  1 

-  721 

30  -  70 

23. 2±  1 

-  721 

co2 

6.6633 

6.645 

77.1 

121.5 

1 . 30±. 3 

1 . 25± . 3 

-211  2 
-211 2 

211.6 

1 . 1 7±. 3 

-211 2 

Sliwinski,  P.  ,  Z.  Physik.  Chem.  ,  Newe  Folge,  63^  263  (1969). 

2  High  Density  Data,  Michels,  A.,  Hamers,  J.,  Physica  _4,  995  (1937). 

3  Stoll,  E.  ,  Ann.  Physik  69,  81  (1922). 
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APPENDIX  F 


CALIBRATIONS 
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THERMOCOUPLE  CORRECTION,  mv 
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TEMPERATURE,  deg  F 


FIG.  28  THERMOCOUPLE  CORRECTION 
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FIG.  29  GAS  CHROMATOGRAPH  CALIBRATION  FOR  n-BUTANE-CARBON  DIOXIDE 
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FIG.  30  GAS  CHROMATOGRAPH  CALIBRATION  FOR  i -BUTANE-CARBON  DIOXIDE 
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FIG.  31  GAS  CHROMATOGRAPH  CALIBRATION  FOR  i -BUTANE-ETHANE 


